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This report presents a prediction method for high-pressure, 
high-euthalpy coretrictor arc heater performance, the conceptual 
desian of S Mw and 40 MW arc hoaters, and supporting documenta- 
tion, An existing computer code for arc heater performance was 
modified by upgrading the radiation model, the thermodynamic and 
transport properties, and the turbulence model. The radiation 
properties model was modified to include visible and infrared 
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atumic bines, the ultraviolet continuum, ultraviolet bands and 
band syaetems, and ulfravinolet atomic Lines, while the radiation 
transport model aas modified fur an absorbing and emitting : 
gaseous mediums, Thermodynamic and transport propervies for air 

soveriag the pressurt range irom 1 to 200 ataospheres and the 

lemperature ranke from 1000°K to 30,000°K were calculated, and a 

turbulence mode) that has been shown to be applicable for devel- 

oping flows, and that gatisties both wall and centerline boundary 

ronditiona, war included, The revised computer code was validated | 

by comparison with existing high-pressure arc heater data, <A 
sealing #tudy uslag the modified computer code was canducted to 
determine the relation between arc fvater performaysce and arc 
feos ken paramgetars, and this information van used in the design 
of 3 WWoand 40 We constrictor arc heaters tor operation in the 
150-200 atmosphere pressure rauge, with maas-average enthalpies 
of 6000-8000 Btu. iba, 
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SECTION } 
IWTRODUCT 10h 


Realistic eiaulation of reentry heating conditione experienced by high 
peiformance reentry vehicles requires the combination of high tent stream en- 
thalpy and high stagnation presaure. Arc heaters offer the potential for 
achieving theses required conditions. The presently available Huele-type arc 
heater provides the required presaure capability but cannot match flight en- 
thaipies. The segmented constrictor arc heater has baeh eployed extensively 
in low-to-moderate pressure, high enthalpy reentry simulation but has not been 
employed at high pressure. Kecent low-co-moderate power teats eat AKC and pre- 
liminary analyses at Aerothers have demonstrated significantly improved enthalpy 
eapability at Kigh pressure fur jhe constrictor arc as compared to the Huele- 
type arc, The heceasary analysis techniques to allow the performance and de- 
sign optiadievation of a high power, high pressure conatrictor arc heater are not 
everiabie, however, Such techniques dre necessary te eliminate ur at least 
minimize the very costly (both financial and achedute) design and hardware jter- 
atlone associated with the eapirical development of auch an are heater. 


Thik report pemsente the develuperent of the me“-.teary accurate analyeia 
touhnique for predicting the resformance atid operating characteriatics of cone 
etrictor arc heaters. Proper physical sodele which are applicable for the cus- 
plete rarge of pressures (to over 200 atm) and other conditiors of interest 
warty incorporated into an existing computer code which was aleo further modi- 
fied for improved capabilities. The resultant computer code was validated 
wie Ugh comparisons of ,  .dicesane with available experimental data. The veli-~ 
dated cade was then employed to determine the relation of performance capablli- 
ties to the various design and operating parameters. Tinally, the conceptual 
design including baait geometric and operating variables was developed for 
woderate and thigh power operation. The performance goal on which the devigne 
vere based wse simultaneous operation in the 150 to 200 atmospheres total 
pressure range and the 6000 to 3000 Btu/lb bulk enthalpy range. 


Tre following briefly describes the report content. Information about 
predictive procedures is d‘a::weed in Sectiona 2 to 5 in terme of previously 
available prediction techniques, iaproved phenomenology modeling for the radiation 
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losses, the thermodynamic and transport properties, and the turbulent model, 
respectively. These ere followed by Section 6 which describes the validation 

of the computer code predictive procedure, and Section 7, which presents the 
results of the acaling study. These ceoctions eummarize essential deteile of 

the technical work, In most cases, additional details are given in a series 

of supporting appendices. Section § presents the conceptual designe for the 

% Of and 40 WW constrictor arc units. Finally, the conclusions of the study are 
presented in Section 9, The aupporting appendices then fellow, che laet being 
a user's manual for the ARCPLO Version 2 code which was developed in part in the 
present study and represents an autumatad version of the predictive procedure. 
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SECTION 2 
PREVIOUS PRODICTLAVE CODES 


Two existing predictive procedures were reviewed for possibie use is the 
present study. The Watson and Pegot (Reference 1: nrocedure was teveicped Jor 
the enalysie of constrictor arcs operating at iow praasuree ant. consequently, 
does not consider phenomenological events which are important at agh pressures, 
as discussed below. This procedure offers the advantage cf seusd coumerics etch 
are suitable for extension to analysis of flows in igh pressure sros. tn addi- 
tion, the procedure offers the advantage of familiacity and has been sant er 
tensively in previous studies to generate predictions which car be aset in the 
present study as baseline data for the evaluation of chanses is the . Sencmpncne 
logical modeling. for instance, the Matson and Peqot code with emripical oor: 
rections has been used by Aerotherm since 1949 for all ar- Seseter desics 
activities. The Cravee and Wells iReference 2} proceture sas the sate ghort- 
comings with recard to aodelina and the same strenqths with retard tc the su. 
merics, but it does not cffer the advantace of hich familiarity to wn ewepet ing 
body of predictions of flows in Sagh pressure arcea. Sased on thence cose, lera- 
tions the Matson and Pegot iReference i: predictive croceduce was astectetd vs 
the baseline for the present study. 

For hich pressure predictions, the chenomenciodicai modeling earicye d by 
Watson and Peqot ‘Reference 1! is inadequate for accurate predics pane of cat ae 
tion flug. [tf also inclodee only iow cressure thermstyc etic ro cryin neee 
perties, and incorporates a someuhat siscliietic curtalest mste whe vata 
properties model does not inciude the folicwias: 

e VYasaible and infrared atotic Lines 

e Ultraviolet continuum 

e Ultravici¢c: basds and basd systems 

e Ultraviolet atomic ines 
while the radiation transpert model eeploys the optictiiy than artroeinarice 


which does not allow self-absorption. The 
diative loss predictions to be lew. while ¢? 
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eredictions to be bigh. At times the teo eppromimations will compensate; how- 
ever, one cannot depend on such good fortune when the radiation flux is the 
dominant loss méchanism (which it ie for high pressure conditions). 


the Wstoon and Peyot thermodynamic and transport properties are subject 
to the following approximations and constraints: 


Air is approximated ae No. dissociated and lonieed 
The pressure iw limited to 1 < p< 10 atm 


es? @e@ 


Thermodynamic proparties are not state-of-the-art 
° Transport properties do not ewpley the moet recent crais sections 
awd their tuxbulent sodel empluys the following idealfirationns + 


@ Amixing length obtainad from the work of Nikuradse (Reference 3) 
ated divided by 2 
@ <A ufaty tucbuvlent Prandel number 


@® An vversimplitied treatment af the effacta of constrictor wall 
roughness 
The thermodynamic and transport property data can be expected to be in substan- 
taal error because they are both aut of date and subject to extrapolation errore. 
Yaw Aurbulent model ig not valid for non-fully developed flows or the flow near 
the centerline of the constrictor tube ar the immediate vicinity of the wall. 


The devetoprent of mare accurate radiation, Cheneodynamic and transport, 
and turbulence modede and data far xanclusfon in an upgreded computor code (s 
presented in the following threr sections. 
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ROCIO 3 
THK RADIANT FLUX IN 9 CONSTRICTOR ARC 


RADIANT THARSPORT 

The accurate calculation of radiation transport within a constrictor 
arc requires considerstion of the gecaetry; namely, # right circular cylinder 
of high length-to-diameter ratio as show in Figure 1. It aleo requires con- 4 
gidetation of the spectral absorption 4nd emission from the gaseous media. 7 
fo obtain the present transport formulation, these features sare combined with 1 
the follewing key aeaumptions: 


@ Wedia ie nonacattering 

® Constrictor walls aro black and maintained at constant tanperature 
@ Cylinder in of infinite length | q 
e Yeaperature coea not vary axially 

e Exponential kernel approximation ia valid 


Theses assumptions restrict the applicability of the analysis to src ficws whic 
do not have appreciable particulate concentrations, which do not have haportant mm 
end wall effects, which have gradiente in the radial direction much larger = 
than those in the asial direction and which do not have walle made of or coated 4 
with reflective materials, None of those ara viewed as serious restrictions 
for the applications envisioned in the present study. 


Consider ea writ aree at Point ¢ (Figure 1) situated at an axial distance 
toand a radial distance rr. Let A-C-N represent a ray having spectral intensity 
1 at Point C directed toward A. For thie system, the epectral flux in the J 
wadial direction is obtained by integrating over a1) the rays passing thiaagh 
C, hie.e 
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Hquation (lt can be vricten ir terag of exponential integral functions 3400 
and Dy ix) (gee Appemlin Ab, which can be approximated by an emponential 
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t Figure 1) Lylindrical geometry end coordinate system. 
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D, (=) © a exp (-bn) (2) 


This approximation allova an analytic integration of Equation (1) ower the # 
variable and results in 


+ - 
a, ir) = a, ir - a,(r) {3} 
whare 
af2 
a, (nr) s Sf com y co’ (r,y) ay (4) 
° 


and where the angular directional fluxes ete, are given in Appendix A. 


Let any of the MN discrete values of the radial coordinate be singled out 
with (he subecript i. he well if located at ivy * Rand the axis of the con- 
atrictor tube at ry = 9. Consider, aa shout in Figure 2, the plane perpendic- 
vler tu the ania of the constrictur tubs, and let f be the index on the redial 
mesh points perpendicular to the axie. Finite difference relations cen be 
obtained (and are given in Appendix A) by assuming logrithmic veriationea Cor 
a with ¢ and for E with +. The anguler directional {lux can be representad 
by @ recursion foimmula which allows significant simplification, i.e., 


. ay , ] 
és bekthed | 
bo lets \ Pili! Fi ae art? ( xi 
tej ) 19,5 i : 
i « ree in — Aad. \ 


oath Baeayy 
where the rn 141.3 ere the optical depth increments. quation (5) includes 
the effect of velt- cabsorption explicitly. 


With known values of spectral absorptian coefficlenta wy gle the 
optical depth incrasante at i-2,3 are generated. Starting at the wall, iem, 


from the known of assumed wail boundary condition, values of Gy are calculated. 


Ove to the symmetry in the ganmetry, ve have oy, 4 s a 1. Avon ine this sym 
metry condition, the ay, j are then computed starting at the amis of the con- 
strictor, se). WAth these calculated quentities, the local spectrai radiative 
Clus qi (ey) may be found fram Bquation ()) cast into proper computational fore 
(Appendin A). 
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To allow assesement of the exponential kernel appromimation, calculated 
radiant flux profiles ace presented in Figure | for @ gray gas ih @ cylindrical 
queometry. The temperature distribution ie esaweed te be iinsar with the radiva, 
and the sbearptica coefficient of the aedium ie assumed to be conetant. The 
calculated radial radiant heat flux distributions are tonmpared with exact cai- 
culations of Reston (Reference 4) and eppronimatye calculations of chibe (Refer- 
wnce 3). Keetun exuployed « numerical integration scheme tc evaluate the ex- 
posential integral functions by (a) and byt), whereas, Chiba and the present 
@ethod used on apyronimation, Chibs (Referance 5) vueed « value of a © 1 and 
be $/4 in the exponential approximation for 0, (kn). whereas, the valuss a © 
Su/l6 and b © 5/4 were selected for the present atudy because they allow addi- 
Rional elaplification in tha analysis, tt ie seen from Figure } thet the re- 
sulte obtained by the present calculational nethod compare excellently with the 
appronimate results of Chiba, and the reeulte are in good agreement with the 
exact calculations of Keaton. 


PADIATIVE PROPERTIES 

The radiative properties of high teaperature air were treated ca a band- 
model basiv. The spectrus wie divided into two gray bands as show in Pigure 4. 
Aboorption qsefficiant data for various preseures (1206 atu) and temperatures 
(1000°R - 30,006°K) were obtained from several sources and are given in Appendix 
a. Roeseland mean opecition were used for the low frequency band, which ie coa- 
sistent with the present, interest in a aelf-abeorbing gas. fhe absorption coef- 
ficients sulected for the high frequeucy bend were selected to correspend to 
the nitrogen ground state photo-lonieation threshold. 


COMPART9O6_WETW MATSON uD PEIOT 

The formulation used in the present study includes the effecte of self- 
abeorption. Thie requires consideration of « specific gecmetry irigkt circular 
cylinder of infinite length) and conmideration of the epectrel nature of tie 
radiation, In contraat, the optically thin appronisation of Wateon and Peqot 
(Reference 1) allows great alaplification in the analysie, in that rediation 
losses need be treated only a6 o heat loss term in the energy equations. Un- 
fortunately, the aptically thin mudel ie not appropriate fur t? 4 high pressure 
constrictor arc enviroement and elould lead to unacceptably high predictions 
(for @ given set of radiation propertice!. 


on the other hand, the Weteun and Peqot (Reference 1) velues of the radi- 
ation properties are not atete-of-the-art and are sleo viewed e6 being incus- 
plete in that all the iaportant cantributions were not included. hie should 
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DISTRIBUTION 
Tlr)-28000~ 20000 (r/R) R 
LEGEND: 

—— KESTON (REF. 4) 
——~ CHIBA (REF. 5) 


RADIANT HEAT FLUX 4, (r), BTU/FT2 SEC 
3 


Cipere } Comper tion of radiant flus profiles, @ «© 0.1 ft. 


Atuc im 6 4GT 


Higwre @. Schematic dlegrem of @ tue-gray-bend absorption coefficient 
sete’. 


causes the prevent predictione to be higher than theire (for a given transport 
formulation) by factora vhich can be ss high es 4. 


the eubstantial differences in both the traneport and properties nodels 
precludes any general statement with regard to which approach gives the higher 
predictions. Indeed, the comperisons which have been made show that the dif- 
ferences can yo either way, Por the impurtent cases of asymptotic flows i con- 
etrictor tubes of 1-1/2 to 2 inches in diameter and in the 100 to 200 atmosphere 
preesure range, the present predictions tend to be about a factor of 2 higher 
than thoee made with the Mateon and Pegot (Refereice 1) method. 
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‘RMTEANIC ABD TRANGPORT PROPERT IEE 


To solve the flow field equations agplicable to a cussetricted arc, the 
following propertion are required in cable format for vee ia the computational 
provedure: 


Therundynemic - +, *%,, y 
Trensport ~ u, &, 9 


in this work, the property tablee aust corer ¢ pressure range of 1 sp 4 268 
ote and a tenpetature range of 1000°« < T < 30,000°R, The priaary weakness of 
the property tables used by Mateon and Pequot (feference 1) if thet thay estend 
wp to only 10 ate. Algo, nitrogen properties ware used to approminaté titted 
Of air, and the nitrogen transport properties used are based upon collision 
Crees-esct ions which are not etate-of-the-art. Due to these deficiaacies, « 
complete updating of the property tables wee a prerequisite to carrying out 
the high-pressure flow fleld analiynes. 


Silewarath, et al. (Aeterence 6) and Gilmore (Meference 7?) provide tabu- 
ler end grephical valuee of thermodynamic properties for sir under the comdi- 
thome of interest. These data vere avkward to vee here becevee of difficaltios 
esecciated with making securate interpolations of the qraphical date for speci ce 
aole frections. In addition, such of these date are presented vith temperature 
and Geagity ae the independent variables; however, presewre and enthalpy are 
the desired independent variables. Therefore, to permit calowlating 7, &, and 
yy for arbitrary valves of the independent veriabies p and h, the celculetional 
procedure [ascribed below ves developed and anployed. 


With regerd to the transport properties, « caiculational procedure vee 
eleo developed even though there are some experimentel data availeble, In pear- 
ticuler, « reasonable amount of experiaente!] date ere aveilable for the eleac- 
teicel end thermal conductivities of a nitrogen plasme (References 6-11), walle 
only limited experimental deta are available for the viscosity (eference 12). 
Yor air, there ace only a few experiaental vaivese of thermal and electrical 
conductivity evailable (References 11, 13), and air viecosity date appear to 
be nonexistent. All of theese data have been ecquired at atmospheric pressure, 
eo that the date can be used to validate ttansport property calculat ional 
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procedures but cannot be used as input toa solve the flow field cquations in the 
200 atm pressure range of interest. A calculational procedure ie therefore 
necessary. : 


A number co’ .snetioc theory calculations have been carried out for both 
nitrogen (References 14-16) and air (References 16-18) plaemes. The calcula- 
tien of altrogen properties by Capitelli and DeVoto (Reference 14] uses the 
best available collision croga-sectiona and ts the aost recent end moat seccu- 
rate; however, only ateospheric pressure was considered. Tha heavily refer- 
unced calculations of air transport properties by You (Reference 16}, Peng and 
Pindroh (Reference 17), and Maneen (Reference 18) have been available for sone 
time, and it now appears that certain collision cross-sections used in these 
treatments are in ¢eriues errs. Furthermore, tho maxiwos preasure conside:ed 
By Yos wae 30 stm, while the other two air property calculations are limited 
AO tempe.,etures below 15,000°RK. For theese reasons, the air transport proper- 
ties were recalculated with the updated modal deacribed below. This model was 
Velidated through extensive comparisons with the one-almoeaphere exper inentel 
date for air and nitrogen plasmas, (t should be noted that all of the experi- 
mental data considered aré very recent, 1970 or later, aid are Viewer) ag being 
the state-of-the-art. 


THERNODYWANIC PROPERTIES 


A chemical mquilibrius computational procedure (the ACE computer program 
(Roferences 1%, 201) wes used to calculate the minture density, enthalpy, and 
apecieos mole fractions for air under the conditions 1 os p ¢ 200 atm and 1LO00*R & 
T < 30,000°R. ‘The ACE code was modified 90 invl de the Debye-Huckel correction. 
(The deteaile of this mxdification sre discussed in Appendia 8.) The Debye-Mickel 
correction is required when ionization is significant to account for the storage 
af potential energy associated with the Coulomb interaction betwee charged 
particles. The net effect of these Coulom interactions je to rafuce the ftoni- 
sation potential, the thermal pressure, and the verious wiature properties in- 
Cluding enthalpy. entropy, density, and internal energy (References 22, 221. 
Under the conditions of interest, the only significant effect is the reduction 
im the itniaation potential, which leads to ahifte in the predicted values of 
charged-particio mole fractions of up to 25 percent. 


The predictions of air thermodynamic properties provided by the modified 
ACH code were cumpared with the values given by Nilsenreth, et al. (heference 
6) and Gilmore (Reference 7?! (ewe Appendix DD. Agreemant on predicted values of 
c and bh was within If percent, wile agreement was alvaye within $ percent for 
the mole fractions of the significant species. 
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‘the new calculations of o and & were also compared with the veluee at 1 
: and 10 atm seed in the Witeon and Pegot procedure (see Appendix 8}. At Comper 
: aturee below 8900R, the Wateon and Pegot {oid) velues of bh are 30-46 percent 

- lower than the new values, while at higher temperatures, they are 14-15 percent 

higher. At the same tuo preveures, there 1s cloee agreement between the old 
and sew values Of density. O72 course, when higher preseures were considered by 
the Wateun and Pegot procedure, the property velues were obtained fram ext rapo- 
lations of the 1 and 10 atm values. It follows that high-pressure propert ies 
determined im this manner can be in substantial error, expecially whan the } 

and 10 atm properties are in error to begin with. 

| 

{ 

| 

j 


TRANSPORT PROPERTIES 
The transport properties were calculated using the eixture rules of fos 
(Reference 16), which are summarized in Appendix 8. Thess etpressions reduce 
to the results of rigorous Ainetic theory in the limit, of a One-apecie yas. 
i Por mixtures, they are approximate in that they exclude the higher order tere 
i in the first Chapman-Enskog epproxiaation (Reference 23). Mownver, calculations 
i based on the simple miuture rules rarely differ from the wore exect first ap- 
prozigation Dy more than a few percent (Reference 16). 


In the Yos formulation, the total thermal Gonductivity K ie the oun of 
tranalational, internal, and ceective contributions. The internal omntcibution 
; ha computed with tho fucken correction (Reference 23), and the reactive thenuil 
conductivity ie based upon the Butler-Brokaw focaulation (Reference 24) for 
multicomponent neutral mixtures which aleo has been shown to be valid for 
pertially-lonized gases in equilibriue (Reference 24}. 


All of the collision integrals (cross-sections) used in the work by Yoe 
wore carefully emamine! and in many cases updated, based on collision integrals 
from References 14, 15, 16. 17. 26. 27, and 29. ‘The detalles of *his invest iga- 
tion are discussed in Appendix ®. for the sake of consistency, the You colli- 
sion integral for # given collicion wee always used when it appeared to be as 
valid a6 that frum eny of the other sources considered, The Yoe collision in- 
tegrais for charge eachange, which make important contrifutions to the reactive 
thermal conductivity. appesred to be too high by «@ factor of up to four. There- 
fore, the charge exchange collision integral for nitrogen was taken fron Capi- 

i Celli and Devoro {Reference 14) and that for oxygen vase taken from Knot, et al, 
(Reference 28). 


tome, 


ray 


Pn eae ey 


Ar ge ag a PE 


seeps a 


AENC TA I$ 47 


Tv) Yor collision integré\s for Coulomb collielone were based on the 
Gvoeduver crove-section aultipiied ty factors ranging froe 6.3 to 12.6, de- 
pending on the particular pair of charged particles. The multiplicative fac- 
tors were obtained through comparison with the el ctrical ani thernal conduc- 
tivities of a fubly-fonized gas predicted by Sritaer and Haire (Reference 29), 
but these latter resulte have been found to be low relative to experimental 
date (Reference 14). Thua, in thie work Coulcabd collision integrals, based 
wpon an unacreaned Coulosb potential with Debye-length cutoff, were taken (inom 
Liborf (Rederencs 27). The Debye length was computed based upon screeming by 
elactrone only, a8 ceocommended by Capitelli end DeVoto (Reference 14). A single 
mltiplicative factor of 0.6, abtained through the compariaona with experimental 
data for electrical conductivity discussed below, wan applied to all Coulomb 
collisions involving an electron. 


The theoretical model for transport properties described abowe wae come 
pared erktically with the avablable experimental data and theories, This 
comparison ia discussed in deteil in Appendix &. Because experimental data for 
the nitrogen plasma are anre extenwive then those for air, the former were used 
as 28 etandard tor comparison, Specifically, the Calculations were compared with 
the ons-atwouphere nitrogen eloctrical conductivity data, and it was found that 
multiplying the collision inteyrala for Coulomb) collisisons involving electrons 
by « factor of 0.6 geve opt iaue agre@eent over the entire temperature range up 
to 24,000°RK. The new sodil chen agreed well with the one-atsospnere resulta of 
Capitelii and UmVoty (Reference 14). Comparisons were alto made with the 100 
atm resulte of Sherman (Reference 16) for To 15,000°K, and gockd agreement wae 
obtained, 


The new model with wodified Coulomb collision integrals was then com 
pared with the available data and theories for the air plasma. Thie comparison 
ie summarised in 7igures 5, 6, and 7. The present calculetions compared with 
the oxperinental deta as woll ae or better than the other available theories 
in ell cases. They are also in good agreement with the one-atecephere resulte 
Of Peng and Pindroh up to fT = 15,000°R, where the latter calculation was ter- 
minated. for electrical comductivity, the present calculations are 26 percent 
hieher than the results of Tos [Reference 16) at temperatures in the vicinity 
of 20.000°K, due to the different Coulomb collistom integrale. The total ther- 
mal conductivity of Yow is up ta 10 percent lower than present calculationa st 
temperatures in the range PCOOR «© FT < 26, 000°R, due to the erroneously high 
charge exchange collision integraia used by Yos. Both the viscosity and the 
total thermal conductivity predicted by Maneen (Reference 18) are in poor 
agreement with the present calculaciones, due most likely to the outdated crose 
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sections used in the former work. Although not shown on the figures, calcula- 
tions were also compared with the 100 atm predictions of Peng and findroh, and 
gond agreement was obtained (see Appendix 8). 


Figures $, 6, and 7? aleo present comparisona between the new properties 
abd those weed by Wateon and Pegot (Reference 1) at both 1 and 16 ata. The 
Matson and Pegot properties are based partially upom the work of Yos for elec- 
trical conductivity and the work of Bansen for viscosity and total thermal con- 
ductivity, and are in substantial error under many conditiona. Consequentiy, 
the new properties should lead to significant iaproveannts in the sccurecy of 
the predictions of flow Clelds within constricted arca. 
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sacTion § 
TURBULENT FLOW MODEL 


For the are conditiona of interest in this study, it is expected thet 
the flow will be turbulent. ‘fhe Reynolds number based on cold flow properties 
and tube diameter is approximately 2 » 10°, a value which far exceeds the usual 
transition value of 3} x 10°, 


By using an eddy viscosity model for turbulent flow, the equations for 
shear stress, 1, and heat flum, q, can be writton as 


t @aotv ¢ a) vv (6) 
wm . oe lak 

q eo « ( 4 -) . (17) 
=. F. dy 


The eddy viecosity, ¢. i& given by 


(8) 


where ¢ is the mizing langth. At the vali, the cining length should satisfy the 
boundary conditions (Referance 30) 


lim le 0 


and (9) 


In the Mateon and Pepot study iReference 1) a modified form of Nikuradee's miz- 
ing length equation (Reference }) was employed (sce Appendix £) which satisfies 
the first boundery condition but gives di/dy + 0.2 aa y + 0. A more suitable 
ecuation for the mixing length in the wall region ia the van Oriest (Reference 
31) Slew of the wall’ model, given by: 
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LAr Le F 
aeaay [ - wo (FE (10) a 
Thies aodel of the mizing leagth satiafius both wall boundary conditions etated 


previcesly, and has been proven effective by other investigators (References 
32, 33). 


All iaforpation on the distribution of the mining length ecross the tube 
readies comes from experimental date (e.q., Reference 34) end evpporte separating 
the flow iato two regione: en lauer region, where a wall nodel for the mining 
length is applicable, end an outer region, where the aixing length ie propor- 
tional to the tube radius. Therefove, the following expreesion tor aixing 
length wes adopted: 


“WTI 70, : 
a “O.4y |) - exe —I.6-5-— for Yo i VSM . 
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where y, i@ a enall distance from the wali and y. i obtained from the continu~ 
ity of &. A comparison of the mixing lengths due to Mikucedee, Watson aad Pegot, 


van Drieat, end the one given in Equation (11) ie show in Figure @ for a typi- 
cal case, 


in addition to changing the mixing length, the turbulent Prandtl number 
used by Matson and Pegat nitede wodification. Matson and Pegot assume a turbulent 
Prandt! number of unity throughout. which ts close to the value often adopted 
for boundary layer calculetians. Rotta (Reference 35) hase proposed the turbu- ' 
lent Prandti number for flow in ducts be given by F 


Pr, 0.05 - 0.08 (E) . aa 


which allows significant deviations from unity near the axie. This value hes 
been used in other recent investigations of duct flows (Reference 36) and was 
adopted in the present study. 


Changing the mixing length model and the turbulent Prandt! number can ; 
have @ significant effect ~w wall heat flun calculations, a6 shown in Figure 9. 
were it be seen that the principal heat transfer mode has been changed from 
rediation to convection for the los -preesure case being studied. At « distance 
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of @5 inches in axial length, the masé-average exthalpy was changed from $200 
Btu/ibm (Wateon and Pegot turbulent model) to 7260 Stu,/lim (present turbulent 
eodel). It should be emphasized that the only differences between the two cal- 
culations presented in Figure § are the mixing length and turbulent Preadtl auwt- 
ber calculations: all other aspects of the two flow field aodele are the came. 


Up to thie poiat, the dlecuseion of turbulence hee aseumed the presence 
of a emooth wall. ia reality, the constrictor wall ie rough, due to both the 
aegnented nature of the wail and the exietence of an cnide scale on the exposed 
eegnent eurfecee (particularly in an air arc}. In both the above fcrmiastion 
end that of Mateos and Pegot, for « emvoth wall the aining length ¢ and, hence, 
the addy viecosity « are sero et the wall. in contrast, for « rough wall Watecn 
and Pegot aseuned the mixing length at the wall was finite and equal to 9.010 
inch for all flow conditions aad constrictor configurations. Nowever, this was 
fekt to be too restrictive in this work. Furthermore, the coughnese associated 
with constrictor segnenta of interest in the present study hae bean seasured 
ot Aerotherm and found to rarely enceed 0.005 inch in equivalent sand grain 
roughness height. 

Yn thie work, wall roughness is nodeled by evaluating Squation (11) above 
at *y ¢ x,° rather than “y*, where K, ie the equivaiont sand grain roughness 
height. This meane that at the wall, y ~ 0, the mixing lemeth ( will be finite 
ant «< 0.4 follows that turbulent camponentse of wal] shear and convec- 
tion meat lus will exlat eince ¢(0) > 6. See Appendix C for further discussion. 


Wall roughness aleo influences the turbulent Frendtl number in the wall 
region. It has been found that réveghness ovgnents wall shear more than it 
augments wall heat transfer, suggesting thst P, in the wail region cen exceed 
unity. tn tale work, Pr, was warited paremetrically and the optiam value ves 
determined to be Py © 3.0 (ene Appendix C). 
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SECTION 6 
CODE VALIDATION 


The improved modele for radiation properties and transport, theraodynamic 
ond treneport properties, and turbulence have been incorporated into the origi- 
eal vereion of the flow field computational procedure developed wy Matecs and 
Peyot (Aeference 1) designated here as ARCFLO Version 1. In addition, further 
minor code aodifications vere perforned to improve the iteration technique used 
to determine the pressure drop for each axial step. ‘The updated code is desig- 
nated here as ARCPLO Version 2. A series of predictions of conatrictor are 
performance vat then mate for operating comlitions wvulere experimental data were 
avéeilable., 


$ ¥ 10m 


The sources of the data available for this purpose are listed in fable 1 
together with the range of constrictor diameter and constrictor lengtha and 
preseures. A listing of all the data is given in Appendix D for the 270 data 
pointe that were collected. 


The following factore wert considered in choneing the best eburee uf ex- 
perieaental data: 

@ igh pressure, high enthalpy levels 

@ Consisteny with other experimental data 

e Sel f-conoistency 


Heximup valves of mass-averege enthalpy and pressure are shvwn in Figure 
10 for the various experiments! data. Lises of constant Hep ere also shown in 
Pigure 10. The AMK conetricted-arc dats is superior to atl of tho other 
eources because it more closely approaches the design goal af 6000-0000 Btu/lbd 
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A power law correlation of all of the experimental data was formuleted 
in order to judge the consistency of the enthalpy and voltege data. Both masa- 
average enthalpy and arc voltege were assumed to vary with current, air mage 
flow rate, pressure, constrictor length and constrictor diameter to some power. 
A miltiple regression computer routine was used to calculate the enponents. 

The following equations were obtained: 
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These equations were used to calculate values of enthslpy and waltage for date 
evelvation. Squationa (13) and (14). while based on extensive data correlations, 
whould be used only for interpolations between given data ranges. They can lead 
to esronedus reeuite when extrapolated beyoud the defining data base. 


Another data test involved the “sonic-flow enthalpy” as calculated by the 
Minsvich iormule (Reference 37): 


tye 2 200 (8), news as 


where A* is the eunic threat ares in equare feet. 


SESYLTS OF DATA EVALUATION 


The date were firet compared with enthalpies and voltages that were cal- 
eulated by neane cf the correlation equations. Of a1) the date, the ASDC qua- 
@tricted arc enthalpy data, shown ia Figure 11, were the best, although excel- 
leat correlations were ateo echieved By the findia dete. The smell amount of 
ecetter i the AZDC data indicates good self-cousistency. 


The voltage comparison of figure 12 shows that the AEDC constricted are 
voltage ie about $0 percent higher tnan predictt? by the correlation formie. 
This discrepancy is epparentiy due to the fact that the electrode witege trope 
are @ larger fraction of the wotel arc voltage for the relatively ehort ARDC arc 
heater. Agein, the small amount of scatter in the voltage data indicates good 
eelf-consistency. 


The final teet of the date from all sourcee ia 3 plot of mase-everage 
enthalpy versus the “soni>-flow enthalpy”. figure 1) shows thie correlation 
for the AZOC enthalpy where the sonic flow enthalpy is calculated veing Squa- 
tion (1S). hen Figure 13 09 compared wilh other such sonic flow enthalpy cor- 
relations, the ABKC constricted arc enthalpy ie superior to all of the others. 


Ae @ reeulit of the abowe data compserisons, runo vere seisected frue the 
ARCC dete and from the Martin Mirietta Corporation dara for the code valida- 
tion. These data, designated as Rune i-16¢, ere listed in fable 2: the final 
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TABLE 2 7 
CXPCKUMENTAL ATA FOR CODD VALIOATION : 
Run hows Vales Sach Tee Vol vet etait piempad 4 
¥ ; \ 821 2080 0.90 17.00 0.055 6,403 76.3 
2 427 RO 0.9 17.90 0. 058 6,0r8 76.0 
3 at 2% 0.9 17.00 0.055 6.909 26.2 
4 ans 3300 0.934 18,00 O.020 5.326 $3.2 q 
Ma $ 370 460 0.95 7.00 Fal 4,$es $1.0 My 
‘ 415 1300 0.943 17.08 6.116 5,96) $3.7 4 
ame iL ? ay} 4730 0.9% 13.00 9.187 4.663 7.6 3 
ry 6) ats 0.9 17,00 0.19 5.200 “4 = 
j ‘ 9 6a? aay 0.9 1.00 0.260 4,448 102.0 fi 
: 10 ear 244 0.96 17.0¢ 0.1% 5 tbh 4.0 
’ 4 629 we 09M 17.00 O.112 $,140 “6.0 
e W $43 ANd 0.99 17,00 0. 12y 5,004 $2.9 
a 635 >obo 9.9¢ 17.00 0,100 6.340 33.9 
i $25 160 v9 17th 0.120 5.02 58.4 
5 454 agno v.9M 17.00 0.283 4.2% 101.5 
16 re) 6176 10" 65.00 o.4a? 10,03? 74.8 4 
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data set selected for code validation were auns 3}, 4, 5, @, 14, and 16. The 
iret five rune include data acquired at the ARGC constricted arc heater faecii- 
ity where preeaures reached 100 atw in a relatively short arc, L/d = 20. The 
Martin Mariett. Corporation data for hun 16 were included to exercise the code'e 
prediction capability for long arce, L/d * 65, where fully-developed or asynp- 
totic flow is obtained. 
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Table 3} eunmarizes the comparisons between experimental data and the 
ARCFLOW Version 2 predictions for both BuIA enthalpy and voltage drop at the 
constrictor exit. The comparisons indicate that the discrepancy between the i 
Version 2 prediction of bulk enthalpy and the corresponding experiental value ‘ 
for a developing are exceeds 10 percent in only one case, while in several 
cases it ie lesa than 5 parcent. Thies agreement is viewsd as being within the 
wuacertainty of the experimental data. The single comparison with a fully de- 
veloped arc ia within 2 percent. The predictions of voltage drop for the ARDC 
test conditions are cunsiatently below the measured values. This is most likely 
due to the fact that the flow field model doee not treat the anode and cathode 
fall regions, far the short ADDS arc, the voltage drops in the electrodes Call 
regions can be 4 significant portion of the total seaeured voltage drop. 


Yor the We arc (Run 16), the wall roughness parameter K, wae parane- 
trically varied from 0.0 toa 0.019 inch, and K, “ 0.0033 inch wes found to pro- 
vide the beat combined prediction of AV and @ when compared to the experiaental 
valves. This value of K, agrees with measurements and estimates made at Aero- 
there. Por the AEOC arc, K, * 9.005 inch was used since the insulator width 
in this arc is somewhat larger than thet for the MC arc. { 


Thy bulk enthalpies are presented in Figure 14 to aliow comparisons be- ' 
tween the Version 1 and Versaion 2 predictions and the axperiental data.® in 
every cate, the Version 7 predictions are superior to the Version } prodictions. 
Consicering only the AEDC data, it if ob#erved ‘hat the Version |! predictions 
are iowers than the measured values, and the dyviations increase vith increasing 
y Preseure, while the tuch maller deviations associated with the Version 2 pre- 
‘ dictions show po particular trend. Murther, the Version 2 predictions for the 
1 jong arc considered in Case 1¢ are in good egreement with experimental data, 
: wubie the Version 1 predictions are substantially too high. In general, the 
i Version 2 predictions compare with the Version ] predictions as tol lows: 


* The Version 2 predictions indicate that a given enthalpy will be 
i reached in a shorter axiai distence 


P ; *The Watson and Pegot version of ARCFLO would not. operate for tun 1S due to an 
, extrapolation of the 1 and 10 ate property tables to negative property values. 
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@ The Version 2 predictions indicate that a lower asymptotic eathalpy 
will be reached. 


A diacucslon of these code comparisons follows. 


Tlowe in short arce are characterized by enthalpy profiles which are 
sharply peeked near the center of the constrictor tube. Energy events in this 
type of flow tend to be dominated by the mixing of the hot core vith the sur- 
rounding cold gases, with turbulent diffusion being the primary transport 
@echansin. Cunsequentiy, the selection for the Version 2 analysis of « turbe~ 
lent Prandtl number which gowee to 0.5 at the center of the constrictor tube has 
the effect of significantly increasing both the predicted transport of energy 
and the predicted axial rate of growth of the bulk enthalpy. 


Man 16 corresponds to « constrictor length for which Cully developed 
or a*ymptotic conditions are approached. In thie particular case, the Version 2 
code calculation predicts twice ss wuch total wall heat flux ae that of Version 1. 
With the wuch lower losses, the Version 1 prediction of ff ia corresponding ly 
higher. Ae diecuesed in section }, the lower prediction of rediative losses by 
ARCTLO Version 1 is due to the fact thet the visible, infrared, and ultraviolet 
lines and the ultraviolet continua are not included in the Watson and Pegot 
model. 


In coucluaion it ia felt that ARCPLO Version 2 provides significantly 
more eccurate predictions in high-pressure applications 39 demonatrated by the 
good agreement between measured values of MH and those predicted by ARCFLO 
Version 2 and the lurge degree of iaprovement relative to the predictions of 
Vereion 1. The remainder of this section is devoted to a brief discussion of 
several physical phenomena predicted by the upgraded version of ARCFLO, 


Figures 15 and 16 present the ARCFLO Version 2 predictions of anial die- 
tributiona for Runs @ and 16, respectively. The anial gradient of @ 19 large at 
the exit of the AED constrictcs, while fo: the auch loager ME conetiicter it 
be nearly aero. This means that higher bulk enthalpies could be sa hieved ina 
the ARK facility if the constrictor vere iengthened and the total voltags drop 
increased while holding mase flow and current constant. In contrast, further 
increases in H in the Me facility cannot be reelised ty elmply lengthening 
the constrictor. 
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Piquree 15 and 16 also indicete that in the thermally developing portion a 
of the flow fleld, the wall heat flux is dominated by radiation. The conves- | 5 
tive heat flux becomes significant only after asyaptotic conditions are ap- 
prosched. fven at this point, convection is typleaity a. wore than 10 pureoat 
ot the totel well heat flux for the elevated operating preseures considered. 
for the ASRC conmetrictor, the wall convection is domineted by the turbulent i 
coatribution due to wall :oughness. In contrast, for the Me case vhere Loth 
bulk Meynolds number and wall rowgkness are smaller, the wall convertion i@ 
apyronimately equally divided between the molecular and turbulent contributions. 
The nature of the radiative and convective wall heat flux predictions in (he ea- 
trance region is e direct result of the entrance profiles considered. The eh- ) | 
trance profiles used in the calculations are discussed below. 4 


Figures 17 and 16 illustrate the radial temperature profiles predicted ; 
by ARCFLO Vereion 2 for Rune 8 and 16, respectively. tn each case, the assumed 
atarting enthalpy profile is essentially the same. the bulk enthalpy corre- 
sponding to the entrance temperature profile is low, being appromimately ¢00 
Otu/lbn. The low we orgy content of the flow at thie point ie assumed to be 
coacentrated in the core; that is, the arc coluen, where significant fonizat iin } 
ie present, resides in a swall region of the center of the flow field, A shore 
distance downstream of the entrance a large temperature spike is generated be- 
cause the Olmic heating is confined to the narrow conducting core of the (how ’ 
field. In both runs, this tCengwrature spike persists past the 17-inch aniel 
position. hen the temperature spike is present, the wall temperature gradient 
is relatively low. As a result, radiation from the core is the major contribu- 
tor to the well hax t flux. Mowever, as indicated for Run i¢, if the flow is 
sllowed to develop, the high-energy core will tend to spread ta the confining 
walle of the conatrictor, and the classical flat profile charactertetic of 
turbulent pipe flow ie approached. Radiation continues to he dominant Ain the 
fully-develuped regime, out the steop wall gradients eleo cause convertion to 
be significent. 
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OCALING STUDY 


The purposes of the scaling study wie to characterize and optimise the nd 
performance of high pressure arc heaters. Specifically, the important parane- 
gare were identified and their effeet on performance eetablished. One of the 
peimary results obtained was a curve relating the maximum mase-average enthalpy 
to preeeure for given values ef maximum pernissible oastrictor wall heat 
transfer rete, Additional oonetrainta such as those imposed by the power sup- 
plies end the test etreem requirements are discussed’ in Section 6. 


The data used for the eculing atudy were obtained fon a series of ARCPLO 
Version 2 computer code calculations. A matrciz of 32 cases vee béentified: 7 
shin metrix is given in Table 4. the input data covers the following renge: al 


Ez 


e 
ae section 7 


Cursent: $00 to 2300 amperes 

Air Mane Plow Rate: 0.125 to 4.0 lim/een _ 
Diameter: 0.75 to 2.939 inchws ’ 
Length: © to 90 inches 


As chown in Table 4, all cases were aucosssfully computed in the firet sttenpt 
encept Case 27. The initial starting eseumptionn for this cease cauped the 
solution to blow up early in the computation aad since the conditions were Rot 
of primary interest a second attempt wae hot made. 
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in order to describe the important trende in the ARCYLO Version 2 per- 
formance data. equations vere sought releting mase-averegs onthalpy, constrictor 
wall heat-treasfer rate, voltage, and efficiency. Theee equations are viewsd 
as ueeful oorrelaticn and iaterpoletion formulae for use in the design optini- 
sation presented in Section @. They should not, however, be vewed to extra- 
polete results beyond date ranges given above. 
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The mase-averege enthalpy wae found to increase with azsial distance at 
@ reletively capid rate to an asymptotic level as show in Pigure 19%. 
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Once the mass-average enthalpy had reached its asymptotic value, furcher i 4 
increases in constrictor length causcl whe radial enthalpy profile to beawoew ; “4 
Platter, DUE did not change the value of the Adee-average enthalpy. Purther, : q 
when the ratio of larval to asymptotic mane-aveyage anthalpy was examingd, 1% 
was found to be primarily a function of the ratio of amial distance to con 
etrictor diambter, %/d, and tolatively Indapenient of conetrictor diameter, 
preseure, alr mage Cliw rate, or current (igure 19). For values of Od greater 
than 15, Che enthalpy-length curve can be apprukinated by 


H g : 
i. fs [> = ann (- sey) (let 
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A eumeaty Vo tle ARCKLO Version 2 values uf mass-average enthalpy, q 
cdestcictor-wall hoat tranafer rate, Soltaye, and current dea given be Table $ 
for a valve of E/d of $1. At this length, HR, | 0.98, 


Correlation equations af the ARCKLOG Version f roaults wore vbtained 
Wain @ MILE Pie eyremsiond SLaLLSLECal Lechiaguae Goa tWase-averaje Galhaipy, 
q ; QUHBECICRO: wahh Neat-transier cate, arc voloage, and efficoncy. The quae 
tions, lor a given vale of afd, ares 
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4 bavetion (21) shows that, for a given pressure and wall heat tranafer 
y rate, the mass-average enthalpy ic solely a function of constrictor diameter. , 
q 4 Curves of maximum enthalpy versus pressure are shown in Figure 20 for several ; q 
= constrictor diameters and an eseuned cometrictor wall heat-traensfer rate of 

a 10,000 Btu/fFt'eec. Thus, it should be possible to attain the “average” 

_ design goal of 7000 Btu/ibs at 175 atmnepheres, providing the constrictor a 
P. dlemeter ies less than one inch. * { 
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1 ree 
: Fractical considerations, as discussed in Section &. limit the general appli- i 
i cation of thie conclusion. For inataace. o 40 1 arc heater should have a 
i constrictor diameter jarger than one inch. 
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SacTtion & 
CONCEPTUAL ARC HEATER OFSIGAS 


The ecaling study discussed in Section 7? provites the basis for develop- 
@ant of conceptual designe for the 5 and 40 MN high pressure, high enthalpy 


constrictor arc heatere, Specific design goals for the arc waite are as 
fartiows: 


e Total mass-average enthalpy: 6900-6000 Btu/lba 
@ Unhanber pressure: 150-205 atm 

@ Minimum operating time: 19 sec 

@ Worzle exit Mech number: 1.7 - 2.3 


A nominal Mach 7 nostle corresponding to an area ratio of 1,79 was chosen for 
design purposes. The maximum levels of current, voltage and input power al- 
lowed for the design are ae follows: 


Parameter 5 te 40_pow 
@ arc current, anps 750 2000 
e arc voltage, kilavolts 1e 30 
@ Input powr, $ «9 


DES LGM _OUIDELIMES 


Thea above performance and operating parameters provide the coviutcainte 
for the designs: there art aleo « number of operating and geometric parameters 
which provide some further depage guidelines. AWaximum values of these ovides 
line paranetera achieved in operational arc heaters serve at least as andicas 
tore bf dvsiyn constraints. The important guideline paremetere are: 


@ Enthalpy-presaurce peramrer, Hep - an indicetor of ovensli arc 
hester performance 


° Conatrictwr wall heat flux, q - an indivator of the cooling require- 
wont s 


6 Are current-caonstrictor diameter paramete:, I/d - an indicator of 
overall loesee and comatrictor heat loed 
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@ Axial voltage gradient, € - an indicator of the maximum constrictor 
disk thickness which is defined by the ailowsble voltege difference 
between adjacent disks, AV 


@ Input power per wait longth, Cl - an indicator of thw local constric- 
tos colum energy ilcadiag 


input power per wait volume, VI/("d'L/4) - af indicator of the over- 
ail constrictor column energy loading 


Coastrictor mess flux, (7) ayy 7 ae indicator of the ccustrictor 
column aerodynamics and ratio of constrictor diameter to throat 
Siometer 


Maximum values of these parameters are presented in Table 6 for the hagh pressure 
experimental deta of the AEOC constricted arc heater and the APFOL Nuela-type arc 
heater, and for all of the data for the ectively cooled arc heatere of Teble 1. 
Consideration of thease results yielded the following maniaun and recommaded 
velues of theaw quideline parameters for the conceptual designe: 


Conceptual Deaign 


cl a etc ct 
® e 


Max iaus 
farameter trom Yable 6 Man ima aecommended 
eB, Beurate “Ibm $2,900 2 . 
Q. Btu/tt' sec 4,620 19,000 5,000 
L/d, emp/ow 638 638 638 
Cc. wolte/ou hs 145 Uns 
av, volts ? 106 100 
Cl, ku/en 210 210 210 
VI/ ind t/4), kwon’ 18.2 Ao 1s 
(ou) gs 1o/ttl eee 16? 250 200 


fven the ainimue performance qos) of 6908 Btu/ibm et 140 atm reguires an in- 
crease of about $6 percent over previously achieved performance. This requires 
in turn an extension of demonstrated capability for sow of the other parameters: 


e Conetrictor wal neat flun, q - roquires high effictency Gboling, 
optimum design constrictor diske 


e Amiel voltage gredient, € - requirow thinner conetrictor dinka toa 
waintain the voltage gradient between adypacent dinke, AV, at ascoep~ 
tabla levels 


pee eee 


“minder deaign goal 74,000 (6000 Btulbm at 150 ste); maxima design 
goal 413,000 (6000 Btu/‘ibe at 200 atm. 
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" TABLE 6 
: OBSERVED OPERATIONAL LOMITS OF WARIOUS AAC HEATERS 
a Ve ko tot eh, 
4 a . ’ ’ 
Ctu-ata T/lem Btustt'sec ampf/om volt/om we bl? cme hu om loft’ sec 
- Conste lctor 
Z | 40,000 4420 %o OS re] ry 1.2 $$ 
‘ : Arron 33,200 00 148 “ Se A) 4.4 167 
mi $2,900 20 638 15 » 210 %$.2 167 
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@® Input power per unit volume, VI/(rd't/s) - requires high efficiency 
cooling, optimum design constrictor disks 


® Conetrictor mass fluz, (50) Lue - emali departure from demonstrated 
acceptable value; results in the requirement for « smaller ratio 
of constrictor diameter to nvrzale throat diameter 


Ic 1Gé SELECTION 


The above guidelines together with the scaling study results of Section 
7 allowed the selection of the optimum conceptual designs. Many conput ations 
were required to develop this optimum design that satisfied the constraints and 
guidelines presented above. Im order to facilitate theee computations, ¢ sin- 
ple computer code which represented the correlation equations for the ARCFLO 
Version 2 results of Section 7? was therefore developed.* fhe results of these 
cosputations, consistent. with the performance goals and operating guidelines, 
were arc heaters with the following basic configurations: 


Arc Heater 


Conf iqurat ian Vartable oe 40 hau 
Conatrictor diameter, im. 0.79 1.75 
Constrictor lengts, in. 2$.5 73.0 
Comatrictor disk thicknesa, in. O.kz 0.20 


Qonatrictor disk spacing 
(ceater-to-center), in. 0.17 0.25 


Note that the design includes a 0.05-inch gap betwee constrictor diske. The 
following paragraphs present predicted performance. 


PREDICTED PERFORMANCE 

The predicted performance af the conceptual designe defined above tes 
presented in Figures 2) through 25 and Tebles @ ani %. The mass-everoge en- 
thalpy es @ function of pressure for both the maximum conditions (q « 16,000 
Btu/tt' sec) and the vecommendad conditions (= $006 Btuvit' sec) as presented 


e 

A lbetiag of the extended BASIC language code utilised is presented in Table ? 
(A .t-ARCFLO). The code applies only to Une results of the ARCFLO Version 2 
code presented in Section 7; it should not be utilized for performance pre- 
ppp peated ontside the range of parameters of the scaling atudy matrix presented 
f Table 4, 
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an Figure <2. iThese resuits aie from the correlation tode which accurately 
characterizes the ASCPLO Version 2% code over the range of conditions of onter- 
est; all other results are directly from the ARCFLO Version 2 code a¢ conde 
tions ted jgecmetricsae close to these of the conceptus) desiun.) From Pigure 21, 
the performance goul can only be achieved at the maximum conditions, and the 

S Wt performance at a given constrictor heat flux level as betrer than that of 
the 40 aw. 


Typical results Crom Uke ARCPLG Version 2 code for a location near the 
downstream end of the constrictor are presented in Tablua 8 and 9, reagectively. 
Those results ara tran ARTKLG Version 2 cangutation Cases JO and @ whach were 
used as oxanples of the radial and axial distribution of progertoes ae pre- 
sented in Pagures 22 through 4. Note that the 5 WM case Cepresetts a more 
fevere Condition than che $0 WW case de.g., qe €200 nmrushe enc ve. $190 Btu,’ 
fy sec. and therefore no conclws ona from quantitative comparikiond are possiade, 


The radial distributions of enthalpy for Both the 5 MM and 40 Mw config~ 
Wratsons ate presented in Fiqure 22, The centerline enthalpy aw almost a fac- 
tor of two higher than the mass-avervage enthalpy, and this facter incrwages 
with uncreasing constrictor: wall beat flux and decreasing caretructor length, 


The axuald dagtrcawutian of radistive and convective constrictor heat (lux 
Ae presented Gn iiquee 23. the radiative flux an by far the dominant tla, the 
convectave flux bebAgd Jesse than $ peroent of the total. 


The constrictor pressure drop and efficiency are presented in Figure 24, 
the effaciency am lower tor the SoM comdituon (Care WO) due to the higher oon- 
atriuctor teat fon and eh lese-than-optamam aac {low rake required by the 
lamated we .tage capedi laity of the abn § FO power aupply. The emeller preaeure 
drep for the SM case ie alec due te the tower flow cate asd therefore lower 


inka flux. 


The net powers ampot - the pamwero te the gan, mo. an presented an 
Fagure 7S ab aw fumctzcn f£ axaal dastance The curve shapo is the same #6 That 
for mass-overage enthalpy Bince concentrated vse aniection at the upstraam end 


of the constrictor weap agbuted for the cotyrutataens. 


A Butmary Of the performance, Iperatano, Qmeometrac, and guideline param: 
ters dor both wre heaters at the recammended conditions ig» $9000 Btu. It! sec 


and 140 atm? we presented twdow: 
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q Design Gosl, Constraint, q 
i &, wtay ito $,900 5,250 $909 to #000 

q Pp, oun 150 150 150 to 200 

4 i, Lom/eee 0.25 4.0 -- , ' 
4 4. Btu/te' sec 5.000 $,000 $000, 10,000 q 
i Vv, ky 9.9 26.3 16 or 30 4 
4 1, anges 690 4,800 740 or 2090 
@ a, in, 0.70 1.75 as : 
' L, in. 25.5 7$..8 -- q 
: €, volte/en 12) 133 119/175 

) Av, volte «“ ry 100/700 , 
4 MiB, Btu-ate “bn 72.300 64,300 -- 4 
: 1/4, smps on ye 8 +56 
a th Reso ”y 199 240/210 

§ Vie Ind be), kw fem! 3 13 As yao : 
(PU) gr Mba te Feme * 240 204,250 i 


q Yor reference, these $ FW and 40 9 comditiona correspond to throat diaueters 
7 of 9,49 and 0.75 inehwon end co exit dianetere of 6.25 and 1.00 Lnctws Coe the 
{ exit Hach number of ¢, respectuvery. Mote that none of the maxiaum guideling 
-, perametc cs presented previously are exceeded. Also, operation at the oondi-~ 
q tame presented an Figure 26 up to flum levels of 16,9000 Beusft sec and 200 q 
ata yialde acceptabihe (but an pom cates maximum) values of the gquidaling { 
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The conclusions derived from the program and the recomendations for 
additional effort are summarized below. 


: @ Accurete characterization of the performance and oparating character- 
: ; istics of constrictor arc heaters, particularly at high pressure, P 
a : requires proper state-of-the-art modeling of radiation, thermodynamic 
ii and transport properties, and turbulence. 


ns: @ Radiation properties must include contributions from continue, 7 

s : linge, and bade for the complete spectrum, and radiation transport ‘ 
7 : muat consider self-absorption; thermodynamic and transport models 

s ~ AMS iaclude proper treatment of charged particles: and the turbulent 

q transport sodel aust adequately characterize physical events, includ- 

log the effects of oomstrictor wall roughness. 


@ Valid apprusimetions and techniques ate available to reduce computa- 

tlowal complenity for radiation without compromise in sccuracy: these 
1 Laciude a tworebsnd rediet»on properties model, exponential approx ima- 
“4 tion of radietion transport, and the use of recursion formulas. 


j e A new computer code, ARCKLO Version 2. which incorporates these pro- fl 
4 per nodele and is based on the procedure of Watson and Pegot (kefer- q 
4 ence 1) Kas been developed and validated for high pressure ias wel) 

{ as low and modereta preseure? conatcictor arc heater applications. 


e This new Code, relative to the original procedure, pradicts thet the 
q bulk enthalpy incresees at a sore rapid rate with aaial dsatance, but 
; reaches @ lower value cf the asymptotic buik enthalpy: the max isun 
) practical constrictor length was found to be defined by a constractorc 
length-to-di ameter catan of 46. 


e Redsation :*+ by far the duminant thermal iose sechenian at high pres- 
4 sure. 


e For the high Reyovlds nanbers typical of high-pressure arce, wall 


q roughness aignificantiy affects wall sheer and heat transfer; further 
: characterization 16 required. 
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APPENDIX & 


% : MOMGRAY, MOMMOMOGANEOUS RADIATIVE TRANSPER : 
x : Im A CONSTRICTOR ARC i 


& The ability to predict the local rediative heat flux is important in the , 
desiqn end operation of a wall-atabilized constrictor arc. Such & prediction i 
; is doubly cvaplicated wwe to the nongray nature of the radiating nedium and due 4 
a to the geceetry. Further complications are encountered vhen the participating l 

; aedium considered is nonhomogensous. Several simplifying aseunpt ions which 
are unresiliatic at high pressure were introduced in the earlier analyses (Raf- 
; erences A-1, A-2). The medium was considered to be: 


@® Optically thin so thet the interlayer absorption mould be neglected 


@ Gray oo that spectral dependency of the radiative properties could 
be ignored. 4 


‘a . In thas analysie the nongray nature of the radiating medium is taken 
Q JNLO account and aleo the radiative proparties ara allowed tu sary spatielly. & 
Moreower, this anetysie is mot lanited eat’: to optically thin or optically 
thick conditions, ‘The loce! rediant heat flux equations are derived from basic 
' primeiples. An exponential kernel approxisation is introduced which siaplifies 
: the radient flux equations, and the resulting equations are then cast in terns 
of an optical depth parameter. A brief description of the numerical acheme is 
given and the results oBtained ave compared with other investigations. 


AMALYS 15 
The governing equation tor radiative transfer in an absorbing and aait- 


ting medium is the equation of transfer, 3.e., 


al 
oe aa ne Oe dandy 


where ® 19 the Planck black body spectral intensity, 1. 18 the spectral in- , 
d tensity traveling along a ray 6 and. is the spectral mase absorption cont fie ; 
cheat correctemd for induced emission. 


Tre spectral radiatave Clux g lft at any radial location ¢ my be eae 
4 pressed ay: 


BAB. 


ft 
ACOC TR 14? 
q H. 
4 a,ir) © fr, cos ¢ da (a-2) : 
: a 4 
¥ where @ is tha angle between the ray and the outward anormal to the cylindrical 
a : eurface and 2 ie the solid eagle. The cylindrical geametry and coordinate sys- q 
x tem are chow in Figure A-1. f 
Bquation (a<l) may be foreslly integrated and substituted iato Equation 
ive’ . (A-2) to yield for a) (References A-3, A-4) 
We . F 
: 72 \ (a'r sing) 4 J cosy 
x aarped / ooe y primes f slyldy ¢ J wly)dy) 
. © © ° 
% ; ‘ 7a 
: at=reinty) 4 cos y q 
‘ * BL yD ly)D, cf vray rs f wlyldy) dy 
a ° ° ° 3 
‘e Wy a 
. cos y r coe y i q 
: + f By lym fy)d, ( { wiy'idy’) dy 5 dy 4 
is ° \ we 
’ 
9/2 \ (Riegointy) 4 qq 
t 74 f coe 1 1B ibd, { f wlyddy) 
q © r coe yY 
fl Vy E 
5 miertainty) y ( ‘ 
\: Wy 
J . f B ty). ty2O, f wly ay?) dy} ay (a-}) ‘ 
J r COs 1 r cos ¥ \ ‘ 
i) 
where : 
} iv 
| ye tet? + plant) (a-4) 
yooeour? - rreanty) 2 (A-$) 
q 
1 ba 
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| } } 
‘ ; 
q 
a AGO TR Na? ® 
ql iq 
aq 4 
‘i \ a 
t ¥ 
bY [° } a . 4 
: Lin) = f axp (- £) a2 (a-6) q 
ig i & L, rr es ( ) 4 
: ; la arriviag at Bauation A-3, it is easumed that the fonscattering sediue a” 
i" : is bounded by @ Wlack eurfece ami is ia local thermodynamic eqilibrius. Fur- of 
a ; ther, At is sncwmed that atial variation of teaperature ie amal) and can be 4 
”" i neyiectod. This approximation 1s consistent with the boundary-layer sinplifi- 4 
‘. ‘ cations adopted ia thie report. ‘ 
5 the Dix) tuactione defined above are known a0 exponential iategrel func~ I 
,* : tiona and are peculiar to the cylindrice! geometry. The Dik) fuactlons have q 
fe the feiluuing properties: 
; $5 O,() © Ohad ed tae?) 7 
tf i and \ 
i. . R* 
q ry q 
» fi “a 
Ai taw 4 
y Dae (xD fraixrax (A~#) E 
q % ‘ 
i" It ke common practice in radiation analyses involving either plane- } 
A paralies geometry or cylindrical weometry to antroduce the exponential kernel ¥ 
4 approximation. aAccoardingsy., folliivwing References A-4 and A-S, we have 4 
4 Dyiny 3 a 97 O™ tAe9) 
y wibenve the comotants 2 and > are swlected auch that they best fit Equation 1A~6) 
q for w= 3. In thie study, numerical values to & and b are assigned to be ; 
: eri 1A-198 
4 and 
r 
7 } 
q be $74 fAodht 
d “5 
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‘ts The local espectrai radiant twat flux 4, (6) ie written as 
ate) © ayy ~ qytr) (aod dd 


¥ where quit) i the radiant {lua ditocted avsy from the location r and q (i) is 
the radiant flux directed tawands the location rf. 


. : The approximate form of the directional epectcal fluxes may be uritton, 
sn terms of anguler directional (lunes G(r,¥), 1.@., 


a 2 
7 Quint = cosy Gi iv.yh ay Kod 3 


: where 

yet 
qq : 
a G*trey > 6a) exp}- (rca - r'sinty) @) © ctr cos yh 

o) 

rhdeenteuntyg 4) 
. Bich exp {+ fee ctr cow HD DE ae 

q a 

a tly cos ¥) 

q ° E(t) exp \- (tle os yt - ei} dt (A=24) 


v 
Gtr.y) © E CR) exp f-(rcam + cteinty) 7) - cir cos wht 
iq toys 1, 4 
rQ mtorr sinty “3 
* git exe }- {kt - tie cos wi} de 


tir cos y) 
(a*2$) 


4 where ciy) 38 the optical depen defined as 


: nly) ° b [over 1A2@) 
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a , Co dyk = 9B tyh (held} 


in thw black body emitaulive power. 


q \ Pauationa (A«12) to (A-27) compiate the formulation of thw spectral - 
i ; : redieat flux equations for a Honhkomoguriequs aedium enclosed in a black-wal lad 
ig comatrictors, it As of jaterest to examine the plysical meanisg of individual 


x : teces in Bavationsa (AvL@) or (A-15). ‘The firat term in Equation (A+14) is the 

La : ‘wall emission that hes deen attenuated by Une 986 medium a8 the cadistion passes 
ae through Points B and C (eee Figure A-1). The second term represents the omis- % 
sion Dy thw yas between Points B and £ attenuated as the radiation pasees trom a 
Be 7 the pulat of emission to Poing C. Radiant energy emitted by the gas volume be- . 
B tween Points E and C, attenuated aa it passes from the point of emission to 4 
: Posmt C ie given by the third term of Equation (A-14). 


7 Analytical eolutions to the above equations are diff icult to ocbtean, ia 
a Reece, 2 MST LGAL belied off devasew, Which £4 SIMPLE, Caputatianaily fast, ( 
and yet accurate. In the following, the numerical method used is dwsecr bad. 


EVALVAT IX OF RADIANT FLUX INTEGRALS ae 


Let the radius of the constrictors be divided anco Nel radnal sumdivasioes. 4 


The wall s6 lowated at rio * Round the axis of the cewstrictor at ron ob. Ae 
thom in Figure A-2, consider the plane perpendicular to tht axia of the cun- -) 
Strictor. het 3 and 1. be the indices on the radial mosh poinge alory the amas A 
and perpendiculer to the axis of Che comsirictor ceapectively. a 


: To eveluste the anguisr directional fluxes So’ trwah, ST ee ut. and Opt acad 
‘ depth - the following procedure :9 adopted. Consider the plane perpemdicalar . 
to the radius verter of any ee As shown in Figura A-d, bee one be tAw angle 

between the radive fant the plane. In evaluaging the optical depen -. toilow- 
tng Niaoalet iMofererce A-6), it 28 aseuteed that the spercgral mass absorpt ps 

q mefticienet . at any alue of y may be written 2° 


AEDE ta S47 


i ; where the quantities y, , andy... 4 are given by 
i Yiog * OL o Eph A {A-iRa) 
- i 
: y ae eee (A-19b) 
( iek.j i+k 3 Se 
q 
i At any value of j, the optical depth increment is, Crom Equationa (A-16) : 
i amd (AV) i 
‘a ‘ - 
' ; 
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Ent in weg : 
¥ ae 
be Combiniay Equations (A-20), (A-24), (A-151, and omploying « loyarit*wic 
‘ interpolation in tems of optical depth for the Slack tody ewiasive power Wie- : 
‘ tribution, the following recutsion rwlatione are obtained for the angular di- q 
% tectional fluxes. At any vaiue of } : 
. a “Ahekad 4 
A wet Z ~ bo 3 Paes e * q 
‘ c* -e@ rytrkey \. aby he Bei i Lerareeaes an paste Cs We sf 
j te) / actd ; we \ f 
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/ Seg eee ON 
; eared en en) q 
taAvdey : 
Starting at the wall. wen. fron sks Anon oundery comdition, values of ‘ 
4 co can be calculated. To evaluate 3° ; the cytaadvycwul symmetry combi tion is 


8.3 
Anvoked. Equetions fasgD) ath (aed 23 ais tw Budetikuked ange Eyuagion (A-b5 


to yteid the tollowing equation tor the directional spectrai fluxes: 
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RADIATIVE PROPERTIES GF GIGN TRMPURATURE ALK 


Rediation properties of high tatpwrature air are complex due to the 
strong vartsation of spectral adsorption vcoefficierat with wavelength ower the " 
Bpectum. The vwarrations in the spectral absorption coefficient are due to 
bound Sree, bound-bouwnd, and free-fr¢ee transitions. 


Detailed calcalatuone of the spectral absorption coetficient are not 
warcanted for this study since they complicate the calculation scheme and also P 
wneveeee The Compuging &.me unvolved considerably. aA sitple band model epprosch 
was adopoed ta chacacterass the varaastion of the abworption cowfficient with 
wave le@ngtd. 


The spectrum (9 to 100 #c! as dovaded atc two gray bands; one band 
covers the range from tlre (Gb ev, and the othor pana extends from 10.5 ev ta 


Sep ev Matha each bend the abeorptaon coed? comme ak, thon, anversent whet 
wavelength. 


Velote of aouvretncn cuedtvcsent fer vavaces presgures and tangeratoces 
are ubtausocu from severe. gevrces, The Aosseland mean free paths fram fohnatan 
kv’ Piates CRelernence ae?) are seed for the low treqoency tend, Covtansam @be 
Borptoos coedforuant wadauys sre atisgned from kheterance [Ae 8) for thw high free 
Poetic tand Yor che Tetperatore cance Poor FOO HR to TOL GOR value cf ate 
korpigon ened ivcaent sve catracted fran enissiuvity gata seportad by Patecman 
ad Mra ak ery ar felezernce A<k:. Figares Asd and A-b show the veriatian of 


ehsctg teat coed tec vert for the tet bands with temperatura fer different pres- 


Soret 
% i] 
ee NC bene model lb be ested sect tne retiative properties are avanle 4 
abe 5 root ere fans ca fatipee. Total cetiatave Jem wt q 


A". Pes ook To wh obtaocte dt be at egret ets bey edtact Mele cue the freqoencsy, 
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Unter the band mote) ass com com Cle lou sd tediative £1ux may Of written as 


q, iri = Yair (A=25) 


where a 36 the total number of bands (in thas study noe 2) and qyir) as the 
radiant flux contributions fram tue [8 nand to the tomal flux whick as given 
by 


q, tps foam a. Lae 263 


' 
ws 


a 


whare he, ag the band-width for the oh? hand. 


hee wy dr) te the lotul tand weagheing fuaceron and detuned by’ 


wolrhs fea fAo27) 


Soe 
wee 


Tiquataon sA-2T) may be ve-owritten an terme of & feactional function wf 
eh faret kird (Reference A-10) onte the mand lauits of the it? bund are Specs~ 
taed, Nete that the local band weighting functson Wir) has numerical values 
tertweernn 0 aod ft. Combining Eqaatsans ied), (a«J6), and (a-1$) and subset nta- 
taog ante Equation tA-2@) leads to the necessary aquatics for the flux fron 
the 1°? band, 


FESULTS AWD U1SCU SSS cm 
Aeduant beet flux dustratalaonk an ow Ceeandricald wediwn ave caic. cated 

for the case cf & fray Gaf with ao singe Sand The Lonmper atece Gust ratwt aor 

28 aS weaMed tr te linear with raduus and the ateworprian ooefficnent cf thaw madi- 

an th Rebate to be anafern. The caldcurated ~adial radiant heat flax dastrae 

tutacne age stow. on Magace Bet. amd are oomeared with exact cadcslatrans .f 

Restor Metevence Ah-} Heston ampli, sc. 2 lame tacd wntegrataau. scheme tc 

CUES de The VeEUNeNta a Vtegrad fame tions Gok. where ag, an the present cal- 

Twi blo hy ROMO, AE HEN tamed mata a ee a BOT OE approaarmaticn web 


Pa ire) Sagere k-t compares the rencite cof sche present soneme with the pegwits 


At ail lg ay, tras teenie tmpemetetat Et 
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at Easton (Reference Av}) and Chile (Reference A-11). Chiba used a value of 
ae i and b © $/4 in the exponential kernel approximation for Bind, whereas, 
in the present study a+ 59/16 and b = 5/4 are assigned. It is sean that the 
results obtained by the present calciiativg) method aGty.e Covorably with 
the approximate reguits of Chibe and the results are in good agreement with 
the exact calculations of Aescon (Reference A-}). 


One of the inherent wraknesses in the present method is thet the pre- 
dicted flux sear the axis of the constrictor is iess accurate. One wey to in- 
crease the accuracy is to have a finer redial mesh near the axis of the con- 
strictor. The advantage of the present computational scheme is that the use of 
recureion relations is sch superior compared ta directly evaluating the radi- 
ant flux equations by, aay. & numerical integration schema. The computational 
algorithm is made simple by 2lininating the integrations required ower the 
anguler and radial coordinetes. 


The strepaqth of tie present approach lies in the fact that it can be 
used te predict redial radiative heat fluxes for all optical conditions of 
interest. This aethod can Le used to determine the effect of “seif-absorpt ian* 
of the cold gas near the wall. The present method can be sasiiy extended to 
inchude @ulti-Sand gases and mixtures of gases as well. 
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APPEMDIX 5 
TMERPICOYNANIC AMD TRANSPORT PROPERTIES 


: Both Versions 1 and 2 of the ARCPUOD code req_ire input of thermodymeaic 
7 and transport properties in tabular format, with pressure and tenperature as 
the independent variables. The property tables are arranged in constant- . 
4 preseure groupes, with each subteble fer s qiven pressure uxtending over a wide ,. 
a range of tenperetures. Version 1 accepts data at only two presaures, ani in mo. 
g : the original woul of Watson and Pegot iReference 8-t} pressures of 1 and 10 ate 
/ were considered, Versicn 2 of the code accepts up lo six constant -pressure 

im tables, amd im thas work pressures of 1, 210, $0, 200, 150, and 200 atm were 

7 Qoneidered, over che temperature range 1000*K «© T < 30,000°R. This appendix 

; discusses in detail the esethods used to generate the property tables for the 
Bix pressures of interest. 


B. THE RMODYMAMIC PROPERTIES ‘9% 


The thatrmadynamis propertiag .. oh, and x are Caiculete?d using the Aeros i 
\ therm Chemical fquilibrium .ACE) computer progran ‘References Bet, Be 3), modie 3 
q 4aed to unchode the Debye-Huckel corraction. Shus subsection presents a brief 
gunmary of tiw ACE formulation for equilibriah gas mixturws, Alec, ancorpora- 
f pon of che GebyesMuckel correctiuns tuto the ACE formulation 49 described. 
ov, Fonaldy, tt vem sdtasg predicesons for .. nm. and x as a fumctiom of p and Tt 4 


exe conparad with val ses available on the literature F 


7 Porat, we untuatafies ATE Rewat int .8 sammharstad. Comeader a gas maox- 
tare comocimet af fo sparcsa@s Ho. pot ee de eee eT, in TALS fyaetem there wiki 
PRLST, Of tie getocal oare, A set of anmbegecutenmt eqealstviun reactions. The 
‘ hanttdne ss ocd UE ket k ae Gh weeda ly Gspesd tm the total sounder of species less 
the sour of elautene ae. Yur Gompuratrcnal purposes. 2 set of species ion the 
ByStan Ls Preseiected and the formation reactions of all other epecres from 


this base set rapresont tne inteperdent get of eqoulibrsum reactiomne: 


AEDC PA Ea? 


where the summation i; over the I base species Ki. 1 © 4, 2, ... , I, and the 
ay are stoichiometric coefficients of the formation reactions. The number of 
Dane species, I, is equal to the nusber of clements in the system. The nunber 
of independent reactions aa then equal to J-0, where ) we Tel, Bet, ... , J- 


Mote that Ios J. 


The most stable (equilibrium) state of this system, if st is maintained 
ot Conetant temperature and pressure, is one for which the Gibos free energy 
of the system is at 4 Ainimum (Reference 8-4). Therefore, associated with the 
qanerai formation cesctiog of Davation ({B-1) is the equilibrium constraint 


# a! ace a 
‘ pan i, (e-2) 


tf the gas mixture is ideal and each species subgas follows the perfect ges 
thermal equation of state, then the partial awolel Gibbs free energy (ctemical 
ponent lal) for speches 7 oan the mixture ke gheen by 


a ot ort inp, (e-3) 


é 


tQwations (8-2) and (®- 3) can be comiined © give an tuprekenom for the 
eyuilebrawm comecant for each independent reactian specafiad Dy Bguataon ced) s 
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The remaining I-] equations required to complete the formulation for closed- 
system gas wixtures are obtained from clement contervation equat sada, ‘ 


The equilibrium formulation gust described is based on the assungt ios a 
that the various molecules are noninteracting except for briaf banary ocean ‘ 
Hy teres which are requiced to establish chemical and thermal eywilibriem, That s 
9 ia, for @ given perticle the tiie between oolkivioos ds mech orwater cna the 
4 Chime tnwolwed Ea collisions. Frog another point of wiew., the partache amter- 

i action potentials ard stahl relative fo these mean Chermal eserpuuo. Whe \ 
q restrictions are sppiicable to « lowedenmsaty gas mixture coma bend nt oa latctinae : 
cally mawtral pertacies (i.a., a tdeab mixture of thermally perfect opaaent. 


Particle interaction porengaal« becom important whwoewrr thes are at rong 

eeouyh to lagivencw tha pargacle ovwe a large portaon af ate dratecteas. 7 ; 
CAN OO wey he Gad MURA U8 at reumeky demo, am whack Gade Atte oan Wie a 
tence Gotweas particles os alwawe my snail chat thee are doe the Sores Qasthat ad i 
a adjacent partaciew. Is con aleo notar of charged partaciet seo oureseon sn the 


M WALKS. MINH Che Casall interaction potamtiad woteeann ta chat b car tac kes a 
LS prvigortusmald to the anveraa od phar SOgurat adn Sub btats bly Se we dt ae 
a baw 4 ech Greater range Kham the repeal wpe pom yl Getta dae aK pe 4 
a tuckes whach typicwlly eacwaw ak tte amverae od RAeue kegepatant Auster ce te 
a Shee Biath om greater powuwsr. Ts Rbk work, Goose aed peetettiad ames. s are ume f 
§ Poreant because charged partacles are present. tte Melye- tucked tracy dee 
eersbed below 18 o8ed to treat shia phancewns, Whar gas domestaeh bre go Mage i. 
that even meutral partigie anturactian pomemtiais antiauance the wae gtate, le a 
second @nd Aayher viarwal cchyvectiaans magt de comna tered Goorm es .it. SE Nira ce 
IMedezsence B-4-. Koweres TOA wit Les Of .ALeTe RT Mere mete meer hho lato 


Cause these vueradi oorvect loons too te slr foc aest ”q 
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Bugm tent To be duro ated Ua pat cam MATA SMe pe Sf TE ge tae tt ban : 
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@e matyoeodgic ecale. st oxae polarized cn oa tacroesoupi. se kle P9C1 AE Te ow 
A HOIROCLAR GL wat Chas preoLarsgatasn. Se PEATE SAT Ewes 2d Ge tates tothe 
expense of the eeectrar bombast exer paws. a7 ther widtse tre. cote aC ar 1 
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The reduction in ronigation eneryy can be deriwed an a purely pmacroa- 
scopic) thermodynamic wanner by extamizing the system Helnholt, Cree anergy 
with respect to lonitatian [Reference ®-2), of imicroscopically? vy sclwany 
Poisson's equation for the outantial on the neighborhood of ay ion sucrouadhed 
by electrons (Reference B-6). in wither ume, aii ie Foun Laat the redection 
is a function of the temperature of the gas and the charged speeches nunber 
dena itien: 


J 
st, wo Jha, e the! typi iT y~ Mn yo (W- 7) 
aay 4! gy e m 7 “ 


taysation (Be 7) be weethen in cogs unuti, and a, » OQ for peeckral atom yp, , k 
for simgly-tomaemd ataw fh. ate. 

The effect of thw Lonktatsan potential lowrung on nmuature canpuagu than 
can be treated by antroducisag 2 COrractaon facton tu the eyaitate ian comet ant 
of Bayaethon (6-4) when written for honating reactnane.  Eepatacn (bed) be Chen 
written as 
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Assuming tha tease species foatve cunpersgaead of the adgutray stoma and she towe 
elwotecn, the sarrectaon factor Saas the forn 
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if it he aeeumed that the lowering of the ion.zation potentials of the yt 
specie in the th and toon) 8 ionization etages hae a negligible influence on 
, their reepective partition functions, then Equation (B-10) can be modified to 
ecrount for the ionisation potentia: lowering by aispiy replacing a2 with 2% - 
art, with ary given by Equation (8-7). shen Bquation (8-10) io generalized to 
the bawe specie formulation on which Bquation (6-4) is structured, by writing 


. 
cae = PCIe SRE tapering soma 


i rel eek 20) . 7 
Co eee) an ' 
| t+ ° TT tnd (@-21) | j 
; ae e : 
' 3 
| with each tera oh the right-hand side of Bqustion (::°11) given by Hquation “ 4 
(8-10) with 5 - ay in place of oe the correction fecter given by Equation : a 
(B-%) falls out. ‘ A 
| the Debye-ttickel corrections to the remaining thertodynamic properties : 4 
ate derived in References 8-# and B-?. fach mixture property ¢ ie atéumed to j 
be @ summation of the unperturbed (uncorrected) value plus a consexibution due ‘s 
| to Coulomd interactions: q 
fe q 
beet ty (B~12) ; 4 
a 
Thus. the internal energy per unit volume is given by 4 
re (med) 
\ 


where 


ee 3 

4 n loa ’ Tt : (fs) . >, ! = va ae: 

Ul ce (cx) re ix ak (Bela 
yw 


Equation (8-24) 18 again obtained by solving Possson'’® equatitn for the peten- 
tral distribution sn the neighborhood of a single charged partocle surrounded 
by # spherically cyemetric cloud of charged particier of cpposite sian iRelor- 


ences Bs6, B-7). The Melwholts free energy 19 given by 
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if 
and, since j 
ee. Fl.n, (e-16) : ; 
one can wrine : 
: 
ar 4 
a a oo (a-17) i. 
Once 7. ie known, the correction to mixture entropy can be obtained: 3 
see + 8, (e-10) a 
1 1 Y. 
6. - a a, - fF} . + : (8-19) 


Also, the pressure correction is given as 
P* Py ¢ dp (B+20) 


where, since 


Rese 1004 (21) 
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One can write 


dtp ey) tl 1 


oe re omenenel r e £ > ei ; 
SS EY Hoe RE 1M eee 


a) 
vince r. ie proportional to v * (nee Bquations (8-17) and (B-14) and note that 
RyP/at ev"'). Pinally, the correction to the minture enthalpy i6 given ee 


| (B+23) 


*e 


alles ia. wes a + Gee heat “Pde yah, fem. é oe cathe ee elem 


AEOC FR 14 OP 


where, since 
fe Seuertp (B~24) 


it fotlows that 
meu, 4 ope fo, (a+2%) 


tn Bqvation (8-20) above, Py is the @o-called “thermal pressure. Since 
dp ie directly proportional to Uy and 0, is a negative quantity, it ful lowe 
that the Coulomb inceractions induce a “*negative’ pressure which serves to make 
the total pleeme p¢ceceure emalier than the thermal preséure. tn wont labora 
tore elesmas, however, this correction is usually quite emall (Roference 8-7). 


Orher alecellaneous relations needed ta incorporate the Sebye-tiickel 
oorrection into the ACE code are the aple fraction defiaition, 


: ry 
xy on (b-26) 


whith requires that Oquation (8-6) be rewritten se 


a 
‘wa } 

, * ) Py (B+27) 
yo) 

The minture equation of state is 

Pp R 

Oo. 7 

val ¥ T (B-28) 


{ ané the conversions from per-unit-volume to per-unit-nase are 
t] 8 
a < (B-29) 


Pinalily. in the equation for siature raective thermal conductivity, 
Oquation (8-36) below, the correction to the enthalpy of species } ie required. 
This correction ie defined in the following manner: 
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Combination of Equations (8-29), (828), (8-25), and (8-14) gives 
(B32) 


The above coulc@b corrections have been incorporated into the ACE code. i 
Predictiane of «, h, and K, from the modified ACE code were then compared with : 
tie calculations in References %-@ and 8-9. In the ACE calculatione, elavan : 
epecieos were coneliteread: 


Se ee 
05+ % o*, ao 
mo, wo’ 
e 


Tables B-) through & 5 preeent o portion of these comparisons. fuble Bel ine 
dicates good agreement. between the unmodified ACE pradictione cad thoes of 
Wileenrath, et al., at 2000°K and 1 atm where the effects of Coulcab interec: 
tions are wseentially sero due to the low degree of ionization, Table -2 in- 
dicates that at 1$,000°K and 1 ata, the Coulomb corrections are small and goad 
agreement with the results of Hilsenrath 16 obtained at thie condition, Finally, 
Table B-} indicates that at 15,000°K and 200 atm, inclusion of the Oebye-miciel 
corrections cen alter the ACK-predicted charged particle number dengitios by as 
much ae 20 ercent, and that these corrections should be included to obtain the 
beet agreement with the results of Hileenreth, et ai. In veneral, the predic- 
tions of ACE with the Coulomb corrections agree with those of Hilsenreth, et 
al... to within 1 percent for ¢ and bh and 5 percent for Rye 
Figures 8+1 and 8-2? present plote of ¢ and bh aa a function of f for the 
ein preseures of interest, a8 viedicted by ACE with Coulomb corrections. Alec 
included are the tabulated values at ) and lu atm and the extrapolated values : 


et 200 arn used by Meteon and Peqot iRefererce 8 1). At teuperatures in the 

vicinity of 4000°R, the Matson and Pegot valves of i are 10-40 percent below } 
the ACB vaivee, white at 16,000°RK - 20,000°R they are 8-15 percent higher 
The Watenn and Pequot values of . at 1 and 19 atm are very close to the ACh 
veluea, but their extra, abctine to 200 ate ta up to 28 percent lower than the 
ACE values. 
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6.2 TRAMBPOST PODPERTIES 


The trensport propertics u, Kk, aad * are calculated using the minture 
mules of You (Reference &-10) and the species mole fractions, apecific heats, 
and enthalpies calculated by the aadified Act code described in dection 8.1. 
The You formation requires numerous collision integrals, and the values 
originally weed by Yos have been updated in thie work through a survey of the 
recent Literature. Ateo, the calculations cerried vut in this vork Rave teen 
conpared extensively with other theories end experimental date available in 
the literature. Thia subsection discusses in detai) the various aspects of - 
the transport properties model developed here, “3 


The exprestions given by You for the transport properties of a partially-~ 
ionized qee mixture are the following (for convenience, use of the subsecripte 
i end j here differs from their use in Section B.t}: 


Ll 
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Ber ia | * x Xs Doers) (B- 34) re , 
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In the above expreseionse, @ i the tots] number df species present (equal to J 
in the nomenclature of Section 8.1). 


The internal thermal conductivity given By Equation (8-39) ie the eo- 
called Recken contribution which eccounte for the tranaport of energy stored 
in the rotational, vibrational and electronic excited states of the various 
species, it id aseumed that the transport of this energy is eseovlated with 
the diffusion process, heave the use of Bquation (9°38) wich q = 1. 


The reactive thermal ootdectivity given by Equation (8-36) esccounts for 
the transport of chesival energy associated vith the diffusion of reecting 
epecion in the mixture, wider the conatraint of chemical equilibrium, an air 
under the conditions of interest, three recombination reactions are the princi~ 
gol contributors to energy tranvport ty diffveiun (Reference B-11), 


@e¢ we ok 
ae hy 
270 0, 


Oquation (8-36) ia based upon the formulation of Sutiler and Brokew (Rsference 
B-121, whieh hes been shown to be valid for embipoler diffusion in e@ pertially- 
lonised gae mixture by Meador and Steton (Referenve 8-13). 


in @quation (B-)6). the summation ower 4 ie 8 summation over all in- 
Gependent reactions in the ainture. Thus, comparing with the subscript con- 
venation veed in G@guation (B-1), the reactions t © 1, 2, .... L iw thie section 
are equivalent to the reactions j = 141, fe2, ... . J in Section 8.1. The 
etoichiomptric coefficients in Bquetion (8-36) are written for ceacticn t in 
the belanced form: 
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which ts equivalent to Bquetion (B-1) eritten in the form 
Yes ", * ‘53 M, « 6 (4h) i a 
ia i if 
where j in Squation (8-41) representa ¢ in Squation (8-40) and v 5° m1. Core a 
teepondind to Pavation (2-40), the heat of reaction per anole of reaction * in ee a 
tquation (8-36) ie given by : a 
an, > %; A, \@-42) ; ‘4 

ie j: : 
where fy includes the Coulomb correction given by Equation (8°31), 4 
tn Equation (B+37), the prime on the summation eign denctes susent ion re a 
. , t a 
over all species except the electron, . i 
In Squation (8-38), the collinion integral wiit+9! has ehe physical Pong 
Significance of an effective droge section, with unite of area, for collisions : aa, 
Setween molecules i and j. The collision integral is given formally by (Refer- : x 
ence 8-3). : 
(p.q) ye (Pog) ’ 
mY . “4, ny 18-4)1 ] 
hj 
where 
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(peal oeE fae tf, a ae te? ? j 4 
(Or 5” rigid sphere ey tt [. oe oc Sad «a3, (B-46) (4 
a. 
2 : a 

v te (8-47) ; 


ae 
re (40) 


end the gaw “inetic cross-section is given by 
ipr p ' ‘ “4 
ay ig) = 2a (t + coos) re {>.q) win a dy (a- 49} ia 


were 2), le the differential crose-section for collisions between aolecules : 
i and 4. _ ie the ecattering aagle ip the center-of-mass syetem, and g is the : 
felative velocity between the coliiding Ooleculea., Gpuation (8-45) specifies 4 
an averaje of the gie-hinetic cross-section weighted by a moment of the Mar- e 
Mallion velocity distribution, In Equation (8-46), @., ie the mean diameter 5 ee 
of eulecules i end § aveuming they are rigid spheres. With the collision in- ¥ 
tegral defined in this manner (Bquetion (B+43)), it reduces to the collision oi 


crose-section ares ‘di, if the tw) particles are actuslly rigid spheres. ee, 


Evaluation of the gee-kinetic cross-section given by tquation (8-49) re» 
quitea knowledge of the intermolecular potential between molecules | and j. Eo 
Gince the evattering engle 1. i8 a function of this parameter (Reference B-$). 
Once the interaolecular potential is known, either from exporinental data or 
@ theoretical model, Equation (8-4)! can be evaluated for the collision ints- 
gral. In the approximate minture cules specified by You, Byvations (@-32) 
through (0-37) in thie work, only the collision integraia for p+ ye | and 
Park are required. Be 


References 8-10, B-11, B-14, B-1%, Bold, Be1l?, and B-10 were consulted 
for collision integrals for the air syetem. Plote of the data for anise! for 
all cullitelone cacept the coulomb collisions revealed thet 
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to within the scatter of the data, where AT, end sf, are comatente. For the 
eake of conbietency, the collivion integral for malécules i and 4 «eed by tos 
(Reference 010) was eleo used here whenever it wae eubetantiated by the veluee 
given by the other references, miuwever, the coblision integrals for charge 
exchange used by foe were found to be too high by @ factor of up to four. 


Thee, tn Chin sdek the Hitteqen chatee onchange integrals ware taken fron 


Capitelli aad Sevoto (Meference 8-14) and thoee for daygen were taken fron 
Manet, et i. (Meferance 8-19). ‘Table B-4 summarizes the constants a}, and sf, 
for #12 bet the Coulomb edllisions. Constants for collisions betweeh @ neu-+ 
tral particle ead & second ion were wot considered, since the nuaber densitics 
for thebe twa eneciea are never sisuitanedualy significant under conditions af 
laveveat, 


The Yos collision integrals for Covloah coliiaions wore based on the 
Gveadever crose-section aultiplied by fattore ranging from 0.) ta 12, 
depending on the particular peir af charged particles. The aultiplicative 
factors were ubtained by You through comparinon with the electrical and thermal 
conductivities of a fully-ionized gas predicted by Spiteer and Hirn (neterence 
B10), But Chese latter reeulte have been found to be low relative to axperi+ 
oentel date (Rafetence b-14). Therefore, in thie work the Coulom collision 
inteqrale were taken from Liboff (teference f-1)), who vaiculated the integrals 
ateuwsing an unecreened Coulumb potential with Oebye-lengeh cuter’. The Libort 
expreasion is (cea units) 


ag(tedl z otf 23 e par [in (2B) - 0.379] (m-92) 


where 


aa, @' 
40 Se {B-$2) 
and the Debye length is 
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(h-39) 


The Debye-length aseuming ecreaning By electrons only is used, ae retomnended by 
Cepitelii and Devn:» (Reference B-i4). The Covlomb colliaion integral for coi- 
lisione involving en electron was corrected using a single @uitiplicative 
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tector, ae outiiaed below. Ali other Coslombh collision integrals, t.e., for 
collisions between various ions, were obtained directly from Equation (B-51) 
with no modifiocetione. 


tutensive cogparicons between the transport property Godel dedcribed 
above and other aodele and experissatel date available in the liteteture were 
carried out. Teble &-S summarises the theoretical caiculations considered, 
and Teble 8-6 summarizes the experimental data considered. Mote that with the 
exoaprion of the Capieelli and DeVoto calculations, alt of the theoretical 
treateente sre felatively Gated. On the other hand, G11 of the experiaeatal 
date sre quite recent. Thie confirns the appropriateness of the transport 
Property eudel updating perforued here. 


‘the primary purpove in carrying out the comparicone between theories . 
and dats wes to validate the property wodel developed in this work. fhe msjor 
portion af the validetion procedure concentrated on comparisons at one atad« 
ephere, dince a1) of the aitperiadntal data and acst of the theoretical calcula- 
tions in the iiteraturo pertain to thie comlition. Nowever, sevwral comperi- 
sone between the present model and other theories were also performed at 100 
ate. 


the following facts were considered in establishing tha validation 
provedurer 


4. Prom @ traneport property point-of-view, an wt, plaema dows not dif. 
fer much trom an air planma (e.9., compare the two caeloulationa 
pertormed by You) 


b&b. There are considerably more exporimental transport property data 


for ™) than there are for air 


oc, There exiates a recent, thorough calculation of 4", plaema t(ranuport 
ersopertiee iCaprtells and Oevatot. 


Goneidering the above cot itraants, it wae decided that the new trenapurt 
property sodel should first be “Cuned® to achieve opt itu agreeaent with the 
theory and experimpntel data for the *, plasma (at one atavephere). then, 
using the aump “*tuned® formulation, the calculations of the new andel were 
compared with the theory end date for §re air plasma (at one atansphere) . 
Finally, tt was ageumed that ell motifications (o the new modal @¢ one atap~ 
sphere are valid alem et the higher pressures of interest, and thie wae con- 
fiend through comparisons betwren the new calculations and the other theories 
at 100 ata. 
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the “*tuning* of the new model was accompliehed by utilizing multiplica- 
tive constants for the various collision integrals. The constants are assuned 
to be independent cf temperature, onepoeition, pressure, etc. This iva fairly 
atendard procedure for forcing agreement between theory and dats for transport 
properties and ie usually required due to the high uncertainty in many Of the 
ecollielon croes-sections, especially théso for Coulomb collisiona vaere the 
shielding proceas is not pieoently well qoantified, in this werk te wes found 
that the only vollision integral correction required was for the Coulieb colli- 
Glens Involving an elect ras. 


Figure ®-) shows the comparisons for the tranepert properties of an Ny 
plasma at one atwoaptere. The frosan thermal comiuctivity Le defined as Yee ) 
Rint (Bquatione (8-34) and (8-35)). The experimental data for elecurical con- : 
Gwitivity were considered to be the primary ctacndard. The calculations of 
Capivelli end deVote were considered to be tho primary theoretiqal standard, 
Mote that Capiteili and deVoto appear to ayree better with the Ny date that 


the other tha@drles coasidured, 


four iterations of the ndw theory were considered: 


a3 


a. Uneoditied cross-sections: without OG ani w’* 


bh. Unmodified cytsa-nections: with O° and N° 
’. ALL Govlosb vollisiot integrals multiplied by 0.6) with Of° and WS 


4. Only Coulomb collision inteyrale involving an electron multiplied 7. 
by 0.6: with o°° and x*?, 


Severe! features of the conpariaons for Ny are evident. 
a. Inclusion of x’° 16 necessary for T < 22,.0C0°R. 


bd. The frozen and total thermal conductivities and the vheactrical con- 
ductivity are quite insensitive ty Coulomb collisions anwoiving 
fonk, eince the thisst ant fourth aterarionm (oe sah do aleve) cive 
enoentiaily the same results, 


ce. The vieccaity is quite inganeitive to Coulomb: collisionn involving 
electrons. for T > 16, 000°R, eince the aocond and third! eter rine 
(hb. and d. above! give essentintly the sant results. 


4. ot followa that o good approach for deterving the multiplicative 
constants 18 to vee the electrical and/or thermal conmdbuctivaty come 
Parison to back out the conatant for electyon-election and electron: 
ton cell.sions, and to use the vitcoeity comperines to back oot the 
eometant for ton-aon collitions. 
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These features also are eatgentially valid for tye air plasma comparisone. 
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Tho finei lteration on the new andel oiovides predictions that agree 
with the w, euperiaental electrical canductivity data to within 10 percent over 
the entire teqperkture range coasidered. in addition, deviations of the pre- 
dictions of the newt aodel frum the ¥, total thermal sondudtivity data never 
axvceed 20 percent for Y <¢ 24,000°K. These particular data @RAIbit large acat- 
ter, and the prediction usually lies within this seatter, Pinally, the now 
eodel predicts M, viscosity within the scatter of the few data points available. 


ror the MW, plagma, the sew model qaneraily compares quits cloeuly with : s 
the rigoroue kinetic theory calculations of Capitelii and DeVoto, being within : Y 
19 percent for total thermal conductivity and ehevtrical conductivity in the 4 
range S000"K «¢ Fa 20, 000°R, and within 20 parcent for teaperatures outside this 
range, The only appreciabla disagreement occure for the vieooaity in the range 
MO, ANG © Pe 10,6004K, where the new acdek prediction ia roughly ¢3 percent 
higher than chat of Gapitelli and Hevota. towever, outside thie tenperature 
range the agteement if better, generally being within 10 percent or less, At- 
teapte (2 reduce tha dladrvepancy for 14,000°K «< fT < 18,000°R were not pursued, 
since experimental data in thie range, which could be used to eubsatantiate 
wither: the new mode] or Capica@li’ and Oewotu, are lacking, 


figure A-4 obows the comparisons for the tranapart progerties of an air 


plasma at one atensphere. The final iteration cf the new aodal provides elec~ , Be 

trical conductivity prédictiang which are within 10 percent af the aaperinenta) 4 

a data for T000°K ¢ 74 1S,000°K and within 20 percent for he unly deta point : 
“ie Outaide this range. The acreement with the total thermal conductivity ia vot - 
i i ae good. being within 30 percent for TOOOTK «© T + 14,000°R and dwviating as a 


woch as TO percent for T © TO00°R. However, 10 thie case there ie only one 
q ‘ met of data with which to compare, and the new model compares with the Jata as 
well aa, or better than, the other theories over the entire teaperature range 


CONS Aer at. 


q in, comparing the thedries for the air plagme, it appears that the new 

4 xxfedl and that uf Peng and Pindroh ore in clone agreenent for all properties ‘ 
for all teeporatures teluw 1$,00078, with the excepticn of the viscosity in the , 
gange 12,0004R © T 5 $4,.000°R. There the new moda! re about 50 percent higher. 
You appeara tu be slightiy low in predicting electrical cecductivity for T © 
22,000°K, Jue to Nis deciei;pe to detersing che marltiphicative constants for the 
Coulomy colligson gntegrala from compariecnm with the predictions of Boatrec 
and HWaen, which are felt to be low themselvae Clapitebi: and Drvatos further, 
fur SODOTR - Po. FHF Poe! peed rot iin of sotal thipmat canductive ty 94 
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Finaliy, Yoo appeare to be substantially too high in hie viscosity prediction 
fer T > 16, 000°, again due to his method of determining the Coulomb miltipli-~ 
cative consteste (thie is also substantiated through the wy conpar isons). 


Te, Baneen prediction for air viecusity ie lower than that of the other 
wodeis for 4000% < T < 10,U00°R, In addition, Hansen's total thermal conduc- 
tivity eppeere to be in grose error for F > 900C°R. 


Figure 8-5 preeents a comparison of the new model with the calculations 
of Sharman for an WN, plaaaa at 10C atm. The agreament between the two viscoe- 
ity calculations jie excellent over the entire tesperature range considered. 

‘The agreement butween the two calculations for frozen and tote] thermal coaduc- 
tivity ia very good for T < 8000°R, tut Sherwan drops below the new model for 
higher teasperatures (although che teaperature-dependent trends are identical). 
Recall) that Sherman's calculation of M, frozen thermal conductivity at 1 ate 
appeata te be low for f > @000°K, relative toe the other theories, including the 
ase aodel and thoes of You and Capitelli and DeVoto. 


Figure 8+é presents « compariaon of the raw godel with the calculations 
of Moneen and Peng and Pindroh for an air plasma at 100 atw. Por viscosity, 
the new model and Pang and Pindroh are within 1) percent for all temperatures 
considered, while Sinsen‘s resulta are generally lower by up to 25 perceit, 
Fur totel tharsal conductivity, the agreement between the new sodel and Peng 
and Pindroh ie excellent, with deviations never exceeding 10 percent. As for 
the ) ate comparisons the Haneen calculation appears again to be groeaely erre- 
neoue. Por electrical conductivity, the new @odel and Pehg and Pindroh differ 
eubstantially for T < @000°K. This is due to the fiict thet the new model uses 
& siqnificantiy lerver e-W, collision integral than that used by Feng and Pin- 
droh. At 8000°R and 100 ate, the acle fraction of M, is 0.48, so that e-N, 
collieions are dominant. At 1 atm and 8$000°K, the aole fraction of ", ie only 
0.06, #0 the eM), collisions are insignificant, thus explaining the good eqree- 
rent between the new model end Peny and Pandroh at thoae conditions. 


Figure 8-7) presente viecumity, Crozen and total thernel conductivity, 
and electricad conductivity for air under the conmditiona 1 s p ¢ 200 atm, 
1O0O*K : T « Fe, 000%R, a9 calculated by the new nodel with corrected elect ron- 
Aon collision integrals. Viscovity av found to be relatively independent of 
pressure for 7 ¢ 12,000°RK. but beccenes increseingly pressure dependent for 
qreater temperatures. VProzren thermal conductivity becomes significantly 
preseure-dependent for T + 8000°K, while a strong pressure-dependence is ex- 
hibited by the total thermel conductivity for tenperetures as low es 3000°8. 
Finally. electrical conductivity 1a a atrong tunction of pressure for elaoet 


Oli temperatures. 
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One noteworthy observation ia that for 12,000°K « T < 13,000°R. abt 
four treasport properties appear to be velatively ineensitive ta pressure 
vartationa for ell properties thie is a ‘cross-over* region below which 
property values decrease with increasing pressure and above which they ine 
crease with increasing pressure. 
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TALE 8-1 


CORPARISOR GF PRESENT CALCULATIGNS WiTw NILSEMRATH, ET &@. 
t + 2600°F, p=) ote 


0. get 0.1984 2 16°" 0.1762 2 10°! 
hy colfgn 47.) at 
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@ @.21 2 16° @.81 a te 
) 0.82 x te" 0.83 » 107? 
0 0.38 5 10°" 0.33 « 19°? 
MOLE FONCTION | C.04 x Yor* Ee 
e 0.29 2 10°" ee 
rat 0.20 2 19°* “ 
0° 0.20 » tor# ee 
wo’ 9.29 « 10°" ee 


TABLE 0-2 


COMPARISON OF PRESENT CALCULATIONS witht WIC SEmMATN, ET a. 
Te 18,080. p+ t ate 


ME with On ACE Uitheat 0-4 


Sarrecttor. _Gereeetion == Mi Npgmreta, of, o). 
D, vee 0.7709 0 1o"* 0.7788 0 10"? 0.7796 3 107° 
: h, caligm 27,44 2? 2 27,409 
a, 0. 9986 2 10°" 6.3706 2 10°* as 
: ‘ : . 
0 0.9018 5 I0°! 0.6208 2 10°! 0.798 2 107° 
MOLE FRACT TOMS a 0.2788 0.2348 0.19 
: e 0. 488 0.9810 0. 
§ a’ 0.2870 0. 2678 0.30 
‘ 0° 0.9701 4 to"! 0.9610 + 10°! 0.61 2 10°! 
no” 0.9097 2 10°° 0.4120 « 10°° oP 
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CONSTANTS FOR EQUATION (8-90) 


(assaminy oily, ta A? dap ttm eK) 


aT) 


+0: 3120 
0.2088 

0. 2676 

6. come 
0.1610 
ff. 2001 
+0, 1963 
0.2500 
0.0860 

0.6798 

0.4768 
6.0" 
9.9978 
0. za 
0.1949 
<0. 1909 
0.2072 
20. 1061 
«, 1008 
0. 2029 
0.2008 
0.1280 

0.6322 
-0. 2 
+0. 2072 
-0. 1009 
0.2072 
0. 2048 
+0. 3979 
-0. 2008 


0.2739 
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1973 
1965 
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1000 - 15,086 
1$.080 
+ 33,006 


TABLE 8-5 
THEORETICAL CALCALATIONS FOR TRANSPORT PROPERTIES ABAILAGLE 36 THE LITERATURE 
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Figure 0-2. Ate density predictions. 
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APPEROIX C 
CALCULATION OF TURBULENT FLOW 


In the calculation of turbulent flows, the ahdar etreas « is composed of 
a iaminar pert and « tuctulent part. By detining an eddy viscosity for turbu« 
lent flow which ie analogous to the kinematic viscosity af laminar fiow, there 
resulte 


y au 
re ote 6 ey ey ' (Col) 
where v * kinematic viscusity 
« «8 eddy viscosity 
» «© #luid density 
tf © mean valocity gradient in the direction normal to the wall 


Similaviy, the heat flux qg is compoeed of laminar end turbulent contri: 
butione ylelding 


wo pty ab 
qe (z . Kay ’ (c-2) 


where ho * thermal conductivity 
voy specific heat at constant pressure 


r. * furbulent Prandt] number 


a * mean enthalpy gradient on the direction normal to the wall 


In the Watson asd Peuut model iReference C+21, the addy viecosity is 
vaicubated by using Prande)  saausga cecath hypothesis, 
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‘ ‘ where it - sisting length. for tlow in amogth pipes, the Mixing length wee found 4 
a é wy Wikuradee (Reference C-1) to be independent of Reynolds nunber for valves of 
he Re > 10°. mikuradse's equation for mixing length ie given in Bquetion (C-é)1 


| feo. - o.00(: - ZY - 0.06 (1 - EJ (e~4) 


i waare 


a : R« pipe radius 
. y “ distance from pipe well 


In correlating data, Satson and Peyot (Reference C-2) found the Mikuredse aixz- 
ing length did not provide good agreement, and reduced it by a fector of two. J 
a Thus, in the Wateon and Pagot eodel, te 1 t.. thie asoumption gave auch f 
= better correlations with low-pressure arc data. " 


a With regards to heat flux calculations, the Mateon and Pegot model as- 
& wuned « turbulent Prandt! number of usity. While thie ie true in the vicinity 
of a wall, it is not true near the center of «a pipe, Mowever,. ho correlation 
prot oma in this regard were noted by Hateca and Peyot. Since recent inveeti- 
gations have found the turbulent Prandtl] eumber deviates considerably from 
unity near the axis for flow in ducte, a turtmient Prandtl number given by 


PL 0.98 = 0.43 (¥) (e=3) 4 


ida wan veed in ARCPLO Veraion 2. 


Mining length formulations which explicitly treat the presence of « 
rough wall do not eppear to be available in the literature. One of the prin- 
cipel ambiguities essociated with this problem is the definition of the ectual 4 
wall location es seen by the flow field when the wall is rough. A second dif- 
ficulty involves determination of the equivalent sand-grain roughness height 
associated with a peculiar roughness geometry (such as segmented conetrictor 
walis), @ necessary step aince most empirical correlations baeed upon experi- 
mental date for wall heat flux and shear sugmentation are expressed in terme of 
equivalent sand-qrain roughnese. 


Order-of -magnitude calculations carried out for the flow/wall conditions 
of interest here indicated that the roughnesse-dominated reuine ia spproeched. 
This meene that the equivalent send-grein roughness height is of the same order 
Of magnitude as the laminar sublayer thickness thet would enist if the well vere 
emocth. for this cease, the friction-factor and velocity-profile deta available 
for low-temperature, inconpressible flow iseo, for instance, Reference C-)), 
can be used to show that the wining length at the wall, t.e.,. et the tops of 
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the coughnees clementa, ia some fraction of the mean roughness element height. 
Uelng thie result tor guidance, it was decided in this work to aodel wall rough- 
nese effects by evaluating the wan Oriest mizing length formale diecuceed in 
Section 5 at “*y + R,” rather than “y*. where y is the distance from the well 
and i, is the equivalert sand-grain roughaces height. At the well, y = o, 

thie then gives ¢) < 0.4 KR, which is coasiatent with the aforementioned low- 
temperature auperimental data base. 


The presence of well roughnese aleo influences thé turbulent Prandti 
manber near the wall. The available experinental date (e.¢., References C+ : ¥ 
end C-4) indicate thet for Reynolds numbers of 10° wall roughness serves to 
eugpent wall ehear by « factor which is up to three times the ourreapondiag 
augmentation of the wall convective heat flux. This is due to the fact that 
the form drag associated with the roughness elements has ano heat conduct ion 
analog. This eleo suguests that Pr, could be an larye ae 3. in addition, the 
detailed profile measurements carried cut in the study described in Reference 
C-4@ involving wall iajection and suction were used to show that tie Rotta cor- 
relation, Squation (C«$) abowe, is quite valid eway from the well. Sowever, 
for y/n < 6.05, #, was found to increase sharply ae the wall vas approached 
and occasionally exceeded even ).0, Saeed upon the calculations described ia 
Bection 6 for ait arcs, in which rough wall effects vere studied paranetrical ly, 
the recommended value for Py, was determined to be 3.0. Por the region 
y/® < 0.08, a linear interpolation between 3.0 ead 0.949, the value given by ‘ 
Equation (C+5) evaluated at y/R © 0.05, was used. . a 
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aAPPENDIz 6 
CONSTRICTOR ARC DATA 


as dincuseed in Section 6, 276 data pointe were gathered fron six 
Gifferent constricter arcs in order to select the nest eppropriate data for 
code validation. A compilation of thie dete is given on tae following pages 
eiomy with material for the identification of each constrictor are facility. 
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Martian Marlette Corporation, Oewver Oivision (0X) 
Gewrer, Colerede 


j fazzle 2 

Corrent, Voltage, SMstricter terest cengtn, Air Flow Mase-Avorege pressure, efficiency, i 

ee oer biamater, teches ates, AAAs ata | percent | 
i L temneen! 

1 380 263 1.0 0.397 = 6.889.001 4,920 1.07 «2.0 

- 2 $3 4160 1.0 0.907 (5.480.188 $,194 8.6 “6 : 

> oO sto? 1.0 0.397 48.08 0.180 8,127 29.9 7.0 : 

4 Oat 1.0 0.997 64.150. 147 10,037 24.06 28.0 

$ 1600 1s 1.0 0.397) 28.33 0,080 14,200 $.14 21.7 3 

6 1200 zy 1.0 0.907 6.500.008 11,208 0.98 24.3 i 

7 400 ase 1.0 0.97 = 13.95 0.084 2,020 0.45 70.4 : 

6 1000 ¥78 1.0 0.997 13.95 = (0.008 15.US 0.085 m4 4 

| ee | 1.0 0.397 «$7.00 (0.082 14,578 0.833 35.1 : 

10688 1.0 0.997 $7.00 0.264 1,97 2.29 $7.4 1 

7) no 4245 1.9 0.997 «= $7.00, 0. Nate &,420 3.63 42.1 4 
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“Total Flew 0.560 lemysecs 0.14) liwvsec. thraugh arc. balance introduced (a plows. 
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APPENDIX & E 
BO USER'S MANUAL FOR ARCFLO, VERSION 2 


Thies appendix provides the information required to operate the ARCFLO 
Verelon 2 computer program. Sections £.1 and £.2 provide input inetruct ions 
aad output descriptions, respectively. section £.) prowides e globe) fiow dia: 
L.. gram and FORTRAN lleting of the code. Section E.4 presenta « sample problem r 
a (the WIC test point discussed in section 6) which was run on a CLC 7600 com- 
puter. Por the sample problem, a listing of the input decke end « few typical E. 
pages of the output are included. .. 


@.1 SMPUT LeRTRLUCT 1ONS 


i tnput to ARCFLO consists of two decks, Deck A and Dech B. Deck 8 con- 
Ae : tains thermodynamic, transport. and radiative property data of air et six dif- : 
ae ferent pressures. Usck & is to be viewed as a permanent deck and no changes 

ti : are to be asde. 

The following ore instructions to assemble Deck A. 4 


a : DECK A (Called from Routine BOUNDC) : 
Card 11 FORMAT (12A6) TITLE 4 


g Title for the perticular cun, used for identification of printed out- * 
i : put. Columns 3-72 are punched with the desired title (alphanumeric). 7q 


Pe ees Bs 


q Card 2: FORMAT (314) RAK, RINC, RTAB 


q ricid 1 (Columns 1-4, RIGHT IUSTIPILO) 
4 BMAX - Naxisum number of ¢xial stations (should not exceed 5000) 


Faeld 2 iColumne $-6, RIGHT JUSTIFIED! 
Kinc ~ Anial station interval for printing output, usually set to 
@ value in the range 30 to 120. 


Fieid ) iColumns $-12, RIGHT JUSTIFIED) 

RTAB - Fleg to print out ingut property tebles (Deck 8) and correspond- 
ing internally-generated tables with finer resolution, leave 
biank for no output, set tc 1 for cutput 

Card 3: FORMAT (214) sA 
Field 1 (Columns 3-4, RIGHT JUSTIFIED) 
MNEEN - Number of radial increments from center to wall, usually set 
to either 13 of 25 


i$) 


. to 
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Cerd 4: 


P| 


Card $: 


Card 6: 


Card 7: 


Card @; 


FORMAT (14) ITURB 


Field 1 (Columne i-4, RIGHT JUSTIFIED) 

ITUMR - Pleg for selecting turbulence nodel, set te 0 for Hateom and 
Pegot model, eet to | for model deecribed in Section 5 of this 
report 

FowuT (14) [START 

Field } (Columns 1-4, RIGHT JUSTIFIED) 

ISTART ~ Flag reserved for restart option (curreatiy sot used, leave 
blank) 

POOMAT (4710.0) AMPS, US, TRCL, B(1) 


Field 1 (Columns 1-10) 
AMPS - Input current in enpe 


Field 2 (Columns 11-20) 
ws - Inlet mass flow rate ia kg/eec 


Field 3 (Columas 21-30) 
TRCL - Transpiration cooling flow rate in kq/eec-a' 


Field 4 (Columns 31-40) 

P(L) ~ Bnlet pressure in ata 

TORMAT (7F 10.0) OLA, THETA, WM, ECRIT, BMAX, RES, TPR 
Field i iColunns 1-10) 

DIA - Siameter cf the constrictor in meters 


Cield 2 (Columns 11-20) 

TUETA - Mogale divergence angie in degrees 
Vield 3} (Crlumns 21-36) 

WH - Wall enthalpy in jowles/hg 


Field 4 (Columns 31-40) 
Ickit — Amiel distance after which current is turned off (i.e@., 
AMPS = 0} in melere 


Vield 5 (Columns 41-50) 

IMAK - Maxisrurm anial distance in meters for which solution is desired 

Faehd © (Columns 51-69) 

82S - Equivalent sand-grain roughness height for constrictor wall in 
meters (0.0000889 m for the MRC arc, 0.000127 = for AEDC erc) 

Field 7 (Columns 01-70) 

The - Turbulent Prandti nvweber at the constrictor wail, generally set 
equal to 3.0 for high-gres¢eure arce 


PORMAT (4P10,0) FIO, EX, EXK, IPE 


isa 


ee 


et 


ARDC TR 164) 


Vaeld 1 iCoilume 1-10) 

F800 — Length of first axial increment divided by the charecteristic 
length £0, ueualily set to FRO - 0.000) (multiplied internally by 
1.@8-06) 

Pisid 2 (Columas 11-26) 

Su - Anial distance increment factor, usually set to ££ = 1.65 

Ligid 3 (Columme 21-36) 

SEE ~ Stability factor, uqually eet to GXZ @ 6.16 


Field 4 (Colwme 31-46) 
GPS — Manion aliowebie relative discrepancy of the mass flow rate, 
apually eet to EPA © 1.0 (multiplied internally by 1.0@-04) 


Card 9: POUT (6710.6) 381, 382, 283. 884, 002, OD) 


These parametere are asecciatead with « code option designed to treet 

vwarleble-area constrictors. This option hes aot heen checked cut and 
gould aot be utilised. Set all 82's equal to SHAK and set all 00's 

equal to OFA. 


Caré (eet) 10; FORMAT (8710.9) H(1,J), J © 1, Meme 


Field 1 (Columns 1-10), Pield 2 (Columms 13-20), etc., eight to a card 
Mi1,3) = Lalet total enthalpy profile in joules/kg (multiplied iater- 
wally by 1.08007) 
Card (eet) 11s POQSMAT (0710.0) O(1,0), Jo 1, MSA 


Field 1 (Columns 1-20), Pield 2 (Columns 11-20), etc., eight to a card 
VI1,3) ~ Inlet axial velocity profile in seters/sec (relative values 
only, corrected to aatiefy glebal mess continuity) 


QBCK ® (Called from Routine wras) 
Permanent deck cards continue. 


.3 OUTPUT OSSCRIPTION 


The AKCFLO Veraica 2 code prints a detailed output block for cach of the 
Cirset three axiel atations. Then, a6 the axiel marching fe continued, eddi- 
tional output blocks are provided at every ximc’” axial station. wote that KIMC 
ie am iaput pareseter. 


Zach output block occupies two pages and coateins both input peranetere 
and quantities which are calculated for the current aziel station. The top of 
the output block conteias the title of the probles which is supplied by the user 
fox identification purpoece. Various input parameters then follow, including 
Ciameter, current, flow rate, wall injection rate, number of radiel anodes, end 
anial etepetze and stability peranetere. The variove calculated quant itics 
eppeer next. Thoee taclude glodal perameters, euch ae bulk enthalpy, and locel 
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parameters, uch as the enthalpy, velocity, and meses flux at each poiat in the 
flow field umere a ncde ie located. In qenerel, the valve of each parameter is 
provided in beth English end 81 unites. 


Yne quantities LOC and 0% shown on the outgut require eume enplanation. 
The quantity LOC is the auaber of pressure iterations required to satisfy the 
tetel masa flow rate at eech axial station. The quantity OW is the errar ia 
the total mass balance, 1.@., 


é -t 
Ow Sale bapyt 
Singur 


where @ is the nase flow rate. 


Towards the bottom of the firnt page of the output block, the current axial 
dietance, masa average enthalpy, wall heat trensfer rates by molecular and turbu- 
lene conduction and radiation, voltage. and efficiency are printed out. 


ta ane version of tre code, a set of diagnostic information is taciuded 
ae the nent to last entry om the first page of the output blowk. The code qu- 
there at Aevothern should be consulted for interpretation o? thie information. 


The final line of output an the firat page contains the input well turbe- 
lent Prandtl! aumber and equivalent sand-grain rovghaees height, and the ceicu- 
lated wall radiation fluxes for the two individual wavelength bands deecribed in 
Section 3, 


The second page of the output block coatains redial distributions of ten- 
perature, TEMPERATURE: mean asbeorption cuefficientsa for the two bande, Kl and Ki) 
amissive power, SEE: Neat Clux potential, PUL (+ S« 47): electrical conductivity, 
RIGMA: gas density, OGNSITY: viscosity, VISCOSITY: mizing length, MIXL: diver- 
gence of the rediative heat flux, OLVQR (2 - 4 - (eqedda Givergeace of the avle- 
vaset conduction heat flum. DIVOC {+ - \ tz (raghs divergence of the turbuler.t 
canduction heat flux, OLVOCT (# - } = ( J); vediel convectica, RADCUN 
tm cy fis amial convection. ANCUN (© cu FED; obese heating, UNMIC HTG (° ont: 


ati redial nasa flux. MMW t= cur. 


a.3 FLOM DIAGRAM ANO COOE LISTING 


Figure E-1 presents the flow dieqram of the ARCTLO Version 2 code. The 
functions of the verious subroutines are briefly described on the flow chert. 
A Portran listing of the code ia presented in Figure E-2. The Fortran veriables 
list 19 given in Reference | and hence is not reproduced here. 
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e.4 SAMPLE PROBLIN 


This sample problem is the solution provided by ARCFLO, Version 2, 

Cor WE Test Point 86> (Num 16 of Bection 6). A Listing of the input cards 
for this problem is shown first. Ilacluded are both Deck A. as described in 
Section E.1, and Oevk &, which is the permanent properties deck for air. Thea 
follows several pages of sample output. Gutput is ikiustrated for the first 
three stations and, in addition, for station 45) {2 = 10.89 inches) and eta- 
tion 16@) (2 © 67.12 inches). ‘The solution waa carried out to station 2500 

(z > 100 i- ches), and the total computer time requirement for a COC 7606 com 
puter was }] seconds (compilation plus emecution). 
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The second sample problem presented is identical to the firse sample 
problem ancept that the distributed mass flow itrarepivation cooling) opt ion 
is utilized, A | vanspiration cooling rate {[TMCL) of 2.00 inate sec 16 assumed. 
ther operating conditions being equal coomparison with the previous cun shows 
the effect of cistributed mass addition on the enthalpy and velonty profiles. 


AC an axial distance of 2.695 inches away from the entrance, che efficiency of 


the arc increases from 0.729 to DO. TH) due to mass addition, and the centus line 
temperature is reduced by about 36$ degrees Kelvin. 
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’ fp fiuid amass flow rate : ; 
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' APPANDIX A 
4 constant used in the exponentiel kernel appromimation, a = */¢ 
4 b conatant used in the exponential kernel approsimation, b * 1.2% 

a Planck blech body spectral intensity 

a cylindricel exponential integral function of order n 

E enissive power 

G enguler direrttonal radiative flux 

i radial index 

r spectral antengity of radiation 

‘ radial index 

' index on the spectral band 

a vocal mu@ber of bands 

s total number of radial nodes 

preasure 

i ape a. 3l radaative flux 

ie Sey as ioral radius 

se compe recor cadiua 

s path lenyth along the line of aight { 

an Loeal band weighting funct son ; 
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the band width 
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7 Stelan-Bbolteman constant 
t optical depen 
bt incremental optical depth 
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al Conetant am Feyeuation 1B-$o7 
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“p, sober specafac heat of specises 4 
4) wean dianeter for hard-sphere mplecule# i and } 
e electronic charge 
P mamture Heimholts free energy, Equation (8-15) 


race 


ated te ear 


‘, oes 5 
: c, partias; silal Gibbs free energy (chemical putentaals of anerien 
2 AN Miah ate 
hie widtie tree eneray cf pure species ) at standard ateate (2) atm 


relative velocity Detween colliding aglecules 


ih elite Senet sia bo eee ae a 
ood 
we 


a siature enthalpy per unit volume 
He heat of reaction per mole of reaction €, Equation (B-42) 
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APPENDIX B (Cont inued) 


: k Bolt rma constant 
L total number of independent reactions ir Minature 
by correction factor for equilibrium constant to account for 


lowering of ionigation potential of specie j, Equation (8-9) 


5. incependent reaction index 

” fixture molecular weight 

*, masa of molecule » 

N represents molecule in Section 8.1) totel num or of species in 
Bection B.2 

ur aumbor denaity of species § in mixture 

Pp mixture total pressure, Equations (B-é) and (B-20) 

Py minature therweal pressure, Fquetions (8-20) and (8-271 

Py part: <1 pressure of species ) in mixture, Equation (B-26) 

Ap pressure correction due to Coulomb interactiona, Equation {B-22) 

a partition Yunct.on for specie } if 2th sonsezation stage 

a universal aaa constant 

s mixture entropy per usit volume 

a mixture entropy per unit mass 

Tt temperature 

UB minature anternal energy per unit volume 
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APPENDIX B (Concluded) 


Vv mixture total volume 

v mixture specific volume per unit mesa 

%, aple fraction of cpecies } in minture 

ay charge number for species j: © for neutral atom, 1 for singly- 


lonized atom, 2 fur doubly-ionized atom, ete. 


Greek 
My constant depending on ratio of aaseea of aolecules i and j}, 
Squetion (8-39) 
a? collision integral parameter, Squation (8-36) 
4 v mixture viscosity, Bquation (8-32); reduced mass, Equation 
a (B-48) 
é wminture moses density 
° winture electrical conductivity, Bquatian (B-+37) 
raped) collision integral for collisions between aplecules i and }, 
Equation (8-43) 
: APFEMDIX C 
; 
ey apecific heat at constant pressure 
fi mean enthalpy 
k thermal conductivity of the fluid : 
‘ winina length 
ty Nikuradse mixing length 
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i {PREMDIK ¢ (oncleded) 
i 
: i. Mateon and reuct mining length » 1/2 te 
3 turbulent Prandt! awaber 
a wall Neat flux 
R constrictor rédive 
a wean velocity ia axial direction 4 
Y distance from constrictor well 
og eddy viscosity 4 
v kinematic viscosity : 
4 
a Fluid density 
t ahear streee 
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